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ABSTRACT 

Aims. We propose a statistical tool to compare the scaling behaviour of turbulence in pairs of molecular cloud maps. Using artificial maps with 
well defined spatial properties, we calibrate the method and test its limitations to ultimately apply it to a set of observed maps. 

Methods. We develop the wavelet-based weighted cross-correlation (WWCC) method to study the relative contribution of structures of different 
sizes and their degree of correlation in two maps as a function of spatial scale, and the mutual displacement of structures in the molecular cloud 
maps. 

Results. We test the WWCC for circular structures having a single prominent scale and fractal structures showing a self-similar behavior 
without prominent scales. Observational noise and a finite map size limit the scales where the cross-correlation coefficients and displacement 
vectors can be reliably measured. For fractal maps containing many structures on all scales, the limitation from the observational noise is 
negligible for signal-to-noise ratios <; 5. We propose an approach for identification of correlated structures in the maps which allows us to 
localize individual correlated structures and recognize their shapes and suggest a recipe for the recovering of enhanced scales in self-similar 
structures. Application of the WWCC to the observed line maps of the giant molecular cloud G 333 allows to add specific scale information to 
the results obtained earlier using the principle component analysis (PCA). It confirms the chemical and excitation similarity of 13 CO and C 18 0 
on all scales, but shows a deviation of HCN at scales of up to 7 pc. This can be interpreted as a chemical transition scale. The largest structures 
also show a systematic offset along the filament, probably due to a large-scale density gradient. 

Conclusions. The WWCC can compare correlated structures in different maps of molecular clouds identifying scales that represent structural 
changes such as chemical and phase transitions and prominent or enhanced dimensions. 

Key words, methods: data analysis - statistical - ISM: clouds - ISM: structure 


1. Introduction 

The interstellar medium has a complex dynamic structure 
on all scales (from sub-parsecs to at least tens of parsecs) 
as a result of various physical processes occurring in the 
multi-scale turbulent cascade of molecular and atomic gas. 
High-resolution and high dynamic range of observations of 
emission line transitions and continuum emission of inter¬ 
stellar clouds p rovide evidence for clumpy structures on all 
scales (see e.g. IStutzki & Guestenlll990l:lRoman-Duval et al.1 
201 11) and anisotropic clouds such as shells and fila - 
ments (e.g. iMen’shchikov et al.ll2010l: iDeharveng et al.1l2010l) . 
Observations of line transitions, total and polarized continuum 
emission provide valuable information about physical condi¬ 
tions (density, temperature, magnetic field) and kinematics of a 
multi-phase gas distributed on different scales. 

The correlation between different structures measured in in¬ 
terstellar clouds, e.g. contours of different chemical tracers or 


the density, temperature and velocity peaks, as a function of 
their size can be used to quantify commonalities and differ¬ 
ences in the formation of these structures. Understanding the 
commonalities and differences seen in different tracers, differ¬ 
ent velocity components or different excitation conditions as a 
function of scale length, will help to infer the underlying phys¬ 
ical processes in the turbulent cascade in interstellar clouds. As 
many processes have characteristic scales, but show their sig¬ 
natures only in particular tracers, it is essential to compare the 
scaling behaviour of various tracers in a same region to identify 
those scales. The characteristic scales could be, e.g., chemical 
transition scales (see, e.g.. [Glover et ali i 2010 ), when compar¬ 
ing molecular line maps of different species, dynamical scales 
for the formation of coherent structures ( Goodman et al.ll998 ). 
IR penetration scales showing up in dust continuum maps at 
different wavelengths (e.g. lAbergel et al.fl996h . ambipolar dif¬ 
fusion scales for the dynamical coupling between ionized and 
neutral particle s, seen when comparing maps o f ionized and 
neutral species (Mc Kee et al.l2QlC)l;lLi et al.l2012h . and dissipa- 
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tion scales when comparing channel maps of individual atomic 
or molecular lines iFalgarone et al.|[l998b . 

To address these issues, we use two different starting points. 
The wavelet analysis has been proven to be a powerful tool 
for detecting structures on different spatial scales in time- 
series and in 2-dimensional maps. It was used e.g. in the A- 
variance analysis measuring the amount of structure in a molec¬ 
ular cloud map as a function of scale and to determine the 
slope of the power spectrum of the cloud scaling (Stutzki et al 


1998). A combination of a wavelet filter ing with the cross 


correla tion function was first pro posed by Nesme- Ribes et aL 
(1995) to study the solar activity. [Frick et all d200ll) improved 


the method introducing the wavelet cross-correlation function 
to study the correlation between galactic images as a function 
of scale. This method and its modifications were successfully 


of astronomical objects and all-sky surveys (e.g., Vielva et al. 

2006; Liu & Zha 

ing 2006; Tabatabaei & Berkhuiisen 2010; 

Roux et al. 2012; 

Tabatabaei et al. 2013). 


significance functions w/( x) and w g (x) characterize e.g. a vari¬ 
able signal-to-noise ratio across the maps or can be used to treat 
irregular map boundaries. 

2.1. Wavelet analysis to measure the scaling in 
turbulent structures 

The A-variance measures the amount of structure in an individ¬ 
ual map as a function of scale, thereby identifying dominant 
structure sizes. In this context “structure” describes the spatial 
variation of measured properties, i.e. a deviation from a flat or 
zero measurement. The amount of structures quantifies the total 
variation in a given map as characterized by the variance. With 
the A-variance, this is evaluated as a function of the size of 
the structures equivalent to the power spectrum (IStutzki et alJ 
1998). This method represents a 2-dimens ional generalisation 
of the Allan-variance method ({Allan! Ill 9961) . The A-variance is 
evaluated from the map / convolved with the wavelet ij/ 


These approaches, however, are not designed for recovering 
the displacement of structures in the data sets on scale-by-scale 
basis and accounting for a variable noise distribution across the 
maps and irregular boundaries. IPatrikeev et al.l ( 2006 ) used an 


F(x , /) = f{x) * i//(x, l) 


by computing its variance 




■x,l)dx\ (1) 


*,y 


anisotropic wavelet transform to isolate spiral features in im¬ 
ages of M 51 for different tracers and analysed their location 
and pitch angle as a function of radius and azimuth but they 
did not cross-correlate the features seen in different tracers. The 
wavelet-based cross-correlation method can be improved by in¬ 
heriting the corresponding formalism from the wavelet-based 
A-variance analysis and combining it with the cross-correlation 
analys is. To account fo r an une ven distribution of noise in the 
maps, lOssenkopf et al.1 d2008al) (hereafter 008) implemented 
a weighting function in the improved A-variance method to 
analyse arbitrary data sets of molecular clouds. The weighting 
function corrects for the contribution of data points with a vary¬ 
ing signal-to-noise ratio as well as allows to perform a proper 
treatment of map edges, independent of their shape, and calcu¬ 
lations in Fourier space and, hence, speeds up the computation 
by making use of the fast Fourier transform. Moreover, 008 
suggested an optimal shape of the wavelet filter for analysing 
the turbulent structures. 

In this paper, we make use of the advantages of weight¬ 
ing functions and develop a weighted wavelet cross-correlation 
(WWCC) method to recover the correlation and displacement 
between structures of molecular clouds as a function of scale, in 
which no assumption about the noise or boundaries are made. 

The paper is organised as follows. Section [2| introduces 
the A-variance and cross-correlation methods. The WWCC 
method is presented in Sect. [3j Application of the WWCC to 
simulated circular structures and fractal structures is described 
in Sect.® Application of the WWCC to observed emission line 
maps of the giant molecular cloud G 333 is presented in Sect.® 
Discussion and conclusions are presented in Sect. ® 

2. Basic theory 

We consider two maps, fix) and g(x), (x = (x,y)), weighted 
at each pixel by a significance function Wfix) and w g ix). The 




=/f ir( *• 


/) - F(l)fdx , 


( 2 ) 


*,y 


where the * symbol represents the convolution and F{1) is the 
average intensity of the map. For uniform, zero average data 
without boundaries, this is equivalent to the wavelet power 
spectrum. 

The wavelet is composed of a positive core, i// c , and a neg¬ 
ative annulus, i// & : 


l//(x, /) = l/f c ix, l) + l/foix, /). 


( 3 ) 


Core and annulus are normalized to an integrated weight of 
unity, so that the integral over the whole wavelet cancels to 
zero. 

For a fast computation of the convolution in Eq. (Q} 
Stutzki et al. (1998) used the multiplication in Fourier space, 
but Bensch et al. ( 200 lb showed that this can lead to consider¬ 
able errors from edge effects due to the incompatibility of the 
implicit assumption of periodic maps in the Fourier transform 
and the boundaries of real observed maps. They suggested a fil¬ 
ter that changes its shape closer to the map boundaries by trun¬ 
cating it beyond the map edges and changing the amplitude of 
core and annulus to retain the wavelet normalization condition. 
While this improves the edge treatment its weakness is that the 
convolution of the map cannot be performed in Fourier space 
any more. 

Ossenkopf et al.l ( 2008a ) suggested a modification that 


overcomes the edge treatment problems, simultaneously deals 
with the effect of observational uncertainties and irregular map 
boundaries, and allows for a computation in the Fourier do¬ 
main. To fix the filter function they increase the map size when 
the filter extends beyond the map edges and apply a zero¬ 
padding to the extended map area. The re-normalization of 
the filter is accomplished by introducing the complementary 
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weighting function w(x) > 0 inside the valid map and w(jc) = 0 
in the zero-padded region. When using a generalized weighting 
function 0 < w(jc) < 1 this can characterize the significance of 
every individual data point, including the effects of noise and 
observational uncertainties. Then the map / p and the weights 
w(x) have to be convolved separately with the positive and neg¬ 
ative filter parts f c (x, /) = /p * 0, M*, 0 = /p * 'I'afa /), 

w c (jc, /) = w(x) * i// c (x, /), and w a (jc, /) = w(x) * if/ a (x, /) and 
the re-normalization is performed when computing the map fil¬ 
tered on scale / 


(4) 

w c (x,l ) w a (x, /) 

When finally computing the A-variance from the convolved 
map, the weights are taken into account in the sum of the vari¬ 
ations (008): 


<n) = 


J[f xy [F(x,/)-F w (l)]W(x,/)dx 


ff xy w F (x, l)dx 


with 


Wf = w c (x, /) w a (x, /), 

and F w (l) as the weighted average of the map 
JJ wp(x,l)F(x , l)dx 


F w (l) = 


JJ Wf( x, l)dx 


(5) 


( 6 ) 


(7) 


The choice of the wavelet filter and its optimisation is 
of great importance for the scaling analysis. For isotropic 
wavelets, the ratio of the core diameter to annulus diameter 
is the critical parameter. 008 found that both, French-hat and 
Mexican-hat filters with a high ratio between the diameter of 
the annulus and the core of the filter, are preferred for measure¬ 
ments of the general power spectral slope, while the Mexican- 
hat filter with low diameter ratios is more suited to sensitively 
detect individual prominent scales. The Mexican-hat filter with 
an annulus/core diameter ratio of v = 1.5 was found to pro¬ 
vide the best compromise, revealing the correct power spec¬ 
trum slope and all spectral features. Therefore, we will use this 
optimal wavelet filter (Eq. (11) in 008) throughout the rest of 
the paper: 




<K(x) 


nl 2 (v 2 - 1 ) 


“ P (< ^ 


\<yU 2) 2 


exp 


-x 

mi 2 


( 8 ) 


where l is the scale of interest. The choice of the Gaussian 
smoothed filter is justified by the simultaneous confinement of 
the filter in ordinary and Fourier space. The usable range of 
scales l falls between 2 pixels, needed to guarantee a reason¬ 
ably isotropic filter shape on a rectangular pixel grid, and about 
half the map size as discussed in Sect. 14.1.21 The method was 
proven to be a powerful tool to characterise the power spec¬ 
trum o f interstellar clouds dStutzki et all 1 9981: lOssenkopf et al 

20089). 


2.2. Cross-correlation 

We can compute the cross-correlation coefficient between two 
images or maps /( x) and g(x) as 


-If 


f(x)g(x) dx. 


(9) 


describing the degree of concordance between the maps. If 
all structures are shifted between the two maps, the cross¬ 
correlation coefficient drops, but the characterization of the 
similarity of the two maps can be recovered when introducing 
a reverse offset vector t while computing the cross-correlation. 
This leads to the definition of the two-dimensional cross¬ 
correlation function as a function of the offset vector, t , as 


C(0 


■// 


f(x)g(x + t) dx. 




( 10 ) 


To account for different absolute scales in the maps one nor¬ 
malizes the cross-correlation function 


C(0 = 


ff x //( x ) - f)(g( x + t)~g) dx, 

<T ./' °~g 


(ID 


where / and g are the means of signals in two maps, respec¬ 
tively, erf and cr g are the standard deviations computed as, 




y/ 2 

JJ(f(x)-f) 2 dx 

, X,y 


( 12 ) 


for erf and equivalently for cr g . The normalized cross correla¬ 
tion function has values between -1 and +1. 

For every given offset vector t , the cross-correlation func¬ 
tion measures the similarity (or correlation) between structures 
in the two maps, when the second map is shifted by t. The 
correlation function for a zero shift, t - 0, recoveres the cross¬ 
correlation coefficient, 


r = C(t = 0) 


(13) 


as the center of the cross-correlation plane. 

The position of the maximum of the cross-correlation func¬ 
tion in the t plane indicates the recovered displacement vector 
between the two maps, i.e. the applied offset for which the two 
structures match best. We compute the displacement vector as 

r = arg max C(t). (14) 

*>y 


The cross-correlation function between two images always in¬ 
tegrates over all scales involved in the map, not distinguishing 
the correlation or mutual displacement of structures at particu¬ 
lar scales. This is addressed by the wavelet cross-correlation. 
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3. Weighted wavelet cross-correlation 


3.1. Theory 

To study the dependence of the correlation coefficient on the 
spatial scale, the two data sets are filtered on scale-by-scale 
basis by means of wavelets and then cross correlated at each 
scale (see e.g jFrick et alJEoOl ). 

Here, we introduce the weighted wavelet cross-correlation 
method to analyse the correlation between maps, / and g, 
weighted at each pixel by Wf and w g , as a function of scale. 
For the filtering, the same formalism is exploited as for the A- 
variance, i.e. Eq. (0 gives the wavelet to filter the maps on 
scale l through Eq. © to obtain the filtered maps F(x , /) and 
G( x, /), and the convolved total weights at each pixel, wp(x, l ) 
and wg(x , /), are computed through Eq. ©. 

To study the cross correlation of the structure in two maps, 
/ and g , as a function of scale, we introduce the weighted 
wavelet cross-correlation function. 


C w (t, l) = 


Cf-g(F l ) 

w,f(J) ^"w,g(0 


(15) 


where Cf-g is the weighted covariance 


Cf-g(T l) 


(w F (x, l) F w (x , /)) x 


If 

x,y 

wg(x + t , /) G w (x + t , /) dx ;, 


(16) 


with 

F w (x,l) = F(x,l)-F w (l), (17) 

and equivalently for G w (x, l) where F W (J) and G w (l) are com¬ 
puted from Eq. (0. The normalization is given by the standard 
deviation 





V /2 

w 2 f (x, l ) Fl(x, l ) dx 


(18) 


of the F w map. cr w G (l) is equivalently computed for G w . The 
normalization WWCC function C w (t,l ) provides values be¬ 
tween -1 and +1. 

The WWCC is therefore a three-dimensional function de¬ 
pending on the offset vector t and the filter size /. In the center 
of the offset plane, t - 0, the WWCC function provides the 
wavelet cross-correlation coefficient r(/), i.e. the degree of cor¬ 
relation of the two data sets on scale /, 


where r x and r y are offsets along v-axis and y-axis, and the 
angular displacement function. 


(p(l) = arctan 



( 22 ) 


For an efficient computation of the WWCC, the convolution 
in Eq. (IT6b can be obtained through a multiplication in Fourier 
space when transforming the weighted functions 


F w {k,l). 


Gw(k, l) 


ff 

If 


w F (x, l) F w (x, l) e lkx dx, 


w G (x, l)G w (x,l) e~ ikx dx. 


(23) 


(24) 


where k = (k x ,k y ) is the wavevector. Using the translational 
properties of the Fourier transform T 


T{w g (x + t, l ) G w ( x + t, /)} = G w (k, l ) x e 


ikt 


(25) 


we can compute the integral for the WWCC function in Eq. 
CE3 through the inverse Fourier transform T 1 {} 


C w (tJ) 


T l {F w (k, l)G w (k, l) x exp ikt } 
^"w,f(0 0~w,g(J) 


(26) 


3.2. Algorithm to compute the WWCC 

The WWCC0 is computed for two images / and g and their 

weighting functions, Wf and w g , through the following steps: 

1. Generate the weighted filtered maps for a given scale l, 
F(x, l) and G(x, /), using the convolution of / and g maps 
with the filter (Eqs. ®, (0 and ®). 

2. Compute the combined wavelet-filtered weights, wp and 
w G (Eq. ©). 

3. Calculate the A-variance (Eq. ©). 

4. Compute the WWCC function (the degree of correlations 
for all possible displacements) on a given scale in Fourier 
domain (Eqs. (1261231241) 1. 

5. Determine the wavelet cross-correlation coefficient and 
wavelet displacement vector (Eqs. 03 and (I2QH221) . re¬ 
spectively). 

6. Repeat the steps 1-5 for every scale, /. 

4. Testing the analytic power and limitations of the 
weighted wavelet cross-correlation 


r(f) = C w (t = 0, /). (19) 

The position of the maximum of the WWCC in the offset plane 
gives the optimum offset vector (for scale /) at which the maps 
are best aligned. We introduce the wavelet displacement vector, 

r{f) = arg max C w (t,l), (20) 

*>y 

which defines the amplitude and direction of displacement be¬ 
tween two maps on scale /. These can be split into the spatial 
displacement function, 

r{l) = (r* + Ty) 1/2 , (21) 


In this section we test the WWCC method by applying it to var¬ 
ious sets of well-defined artificial maps to answer three ques¬ 
tions: 

- Can we determine individual characteristic scales in the two 
maps? 

- Can we reproduce a given displacement of structures at par¬ 
ticular scales? 

- Can we trace a systematic change of the scale sizes between 
the two maps? 

1 The IDL code of the WWCC method is publicly available at 
http://hera.phi.uni-koeln.de/~ ossk/ftpspace/wwcc/ 
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To address these questions, we use two different basic sets 
of simulated noisy maps. In one case, we assume Gaussian in¬ 
tensity profiles of a given size, representing individual smooth 
clumps. Then we have only one well-known prominent scale 
in the two maps to be compared and we can test the detection 
of a possible scale difference and a displacement of the struc¬ 
tures. The other set of data is given by self-similar fractional 


Brownian motion structures (fBm, IPeitgen & Saupel Il988l) . 


They do not have any dominant scale, but are characterized by 
a power-law distribution of scales. They can be described by 
a power law power spectrum of fractal structures and random 
phases in Fourier space. Here, we can test whether any scale 
dependent modifications (smoothing, displacement,...) are re¬ 
covered by the WWCC. The two test cases (circular and fBm 
structures) should eventually represent limiting cases for the 
actual observations of molecular clouds. In one extreme view, 
maps are described as a collection of “spherical blobs”, in the 
other extreme they are described as scale-free, fully self-similar 
fractals. In practice, they typically contain both aspects, i.e. we 
find both self-similar structures and individual characteristic re¬ 
solved structures. 


For these test cases, we investigate the CC coefficient and 
displacement vector when changing the properties of the arti¬ 
ficial clouds, such as the size of the individual structures, the 
spectral index of the self-similar scaling law, and the assumed 
noise level of the observations to test their sensitivity to those 
properties. The map size of all test data sets is 128x128 pixels 
(abbreviated as pix throughout the rest of the paper). 


4.1. Structures having Gaussian intensity profiles 

We test the WWCC method for circular structures having 
Gaussian intensity profiles and some spatial displacement to 
study the effect of the size of circular structures and the noise 
level on the correlation coefficient and the measured displace¬ 
ment vector as a function of scale. 

The Gaussian intensity profile is generated with an ampli¬ 
tude equal to unity and width given by standard deviation cr. 
A weight (significance value) of unity is assigned to each data 
point (pixel) in the map. 



Fig. 1. (a) Maps of two circular Gaussian structures (cr i = 10 
pix and cr 2 = 5 pix) shifted by 20 pix along y-axis. (b) Observed 
with a simulated S/N = 5. 



Fig. 2. Maps of the two Gaussians with S/N = 5 from Fig. [I}) 
filtered on scale of 5 pix (a) and of 45 pix (b). To allow for a fil¬ 
tering in Fourier space, the maps are extended beyond the orig¬ 
inal boundaries before filtering. The extension is treated with 
wp — 0 in the computation. 


4.1.1. Map filtering 


In Fig. [I] we show two pairs of test maps given by circular 
Gaussian structures (with cr 1 = 10 pix and cr 2 = 5 pix and a dis¬ 
placement of the peak by T; = 20 pix along y-axis in the second 
map). As the cross-correlation (Eqs. (ITOl) and (fl5l ) ) only con¬ 
siders the relative offset of structures between two maps, t, the 
result is invariant with respect to an exchange of negative off¬ 
sets in map 1 by positive offsets in map 2. The WWCC always 
characterizes the mutual displacement of structures within the 
two considered data sets. Here, we always displace the struc¬ 
ture in the second map relative to the first one. The upper panels 
show the original structure representing a S/N = oo, the lower 
panels are superimposed by white noise simulating typical ob¬ 
servations with S/N = 5, i.e. <x no ise = 0.2. 


The first step in the analysis is the wavelet convolution (Eq. 0) 
filtering the maps on a given scale. The usable range of the fil¬ 
ter sizes is inherited from the A-variance analysis discussed in 
Sect. 12.11 It starts at a resolution of two pixels. Fig. [2] shows 
the result, F(x, l ) and G(x, l ) for the S/N = 5 case for two se¬ 
lected wavelet filter sizes l = 5 pix and l = 45 pix. The filtered 
maps are larger than the original ones to allow for a convolu¬ 
tion with non-truncated wavelets (see 008). The extension is 
treated with zero weighting in the further computation. We see 
that the maps filtered on small scales (5 pix) are dominated by 
noise, while the structures filtered on large scales (45 pix) are 
basically free of noise. Statistically, the noise will dominate the 
cross correlation for all scales l below a critical size. This size 
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Fig. 3. A-variance spectra of the two Gaussians with S/N = oo 
from Fig. Q](a). The A-variance spectra of the large (left) and 
small (right) circular structures (Fig. |T](a)) are denoted by full 
and dotted lines respectively. Vertical thick lines are used to 
denote the measured prominent scales of the large and small 
circular structures, respectively. 


Fig. 4. A-variance spectra of the two Gaussians with S/N = 5 
from Fig. Q](b). The A-variance spectra of the large (left) and 
small (right) circular structures (Fig. Q](b)) are denoted by full 
and dotted lines respectively. Vertical thick lines denote the 
measured prominent scales of the large and small circular struc¬ 
tures, respectively. 


depends on signal to noise level and will be discussed later in 
this section. 


4.1.2. The A-variance spectrum 


Computing the variance of the wavelet-filtered maps for all fil¬ 
ter sizes leads to the A-variance spectra (<x^(/); Eq. ©). Fig. a 
shows the resulting spectra for the original maps without noise 
and Fig. [4] the spectra for the case with S/N = 5. The A- 
variance characterises the amount of structure distributed on 
different spatial scales. The peaks of A-variance spectra indi¬ 
cate the dominant structure size. For structures without addi¬ 
tional small scale contributions, the A-variance falls off like / 4 
below the dominant size. This is well seen in Fig. [3] for scales 
below about 8 pix. The increasing statistical uncertainty of the 
A-variance for larger scales due to the lower number of statisti¬ 
cally independent entities of large scales in the maps practically 
limits the use of the A-variance analysis to prominent sizes of 
less than half of the map size (008). The huge formal error bars 
in Fig. [3] folio wing 008 are irrelevant here, as they characterize 
the counting statistics of individual structures, being very bad 
for a single structure in an otherwise empty map, while every 
observed map typically contains much more structure. 

The A-variance detects the prominent scales of the circu¬ 
lar structures at Z p i = 41 pix and / p 2 = 18 pix (vertical lines 
in Fig. [3]). This corresponds approximately to the diameter of 
the circles visible in colors in Fig. Q](a), i.e. at a level above 
10 % of the peak intensity. For isotropic structures, the promi¬ 
nent scale approximately identifies the visible size of the struc¬ 
ture (Mac Low & Ossenkopf 1200(1) . We find a fixed relation 
between the prominent scale and the standard deviation of a 
circular structure / p « 4<x. 

When we compute the A-variance spectra of the same struc¬ 
tures simulated with S/N = 5 (Fig.® we find a superimposed 
structure on scales less than 8 pix which is due completely to 


noise. The noise contribution decreases with <x 2 (/) oc / 2 (see 


Bensch et = aD 2001) up to « 8 pix and becomes negligible to¬ 


wards larger scales as discussed for the map filtering above, 
leading to a A-variance minimum at a scale of about 8 pix here. 
At scales above the minimum the spectra behave identical to 
Fig.[3]and reach the maxima at scales of 41 pix and 18 pix. 


4.1.3. The WWCC function 

The WWCC map is finally computed through Eq. (l26l) provid¬ 
ing the correlation for all possible displacements on a given 
scale. It is then used to compute the CC coefficient and dis¬ 
placement vector as a function of scale. The center of the 
WWCC map represents the correlation of the two structures not 
correcting for any displacement. It gives the cross-correlation 
coefficient for the given scale (Eq. (fl9l)). The dominant offset 
between the two structures in the maps is measured as the dis¬ 
placement vector between the center and the peak position of 
the WWCC function (Eq. ®). 

This is demonstrated in Fig.[5]for the pairs of maps wavelet- 
filtered on scales of 5 pix and 45 pix from Fig. [2] The black 
crosses indicate the map center, i.e. a zero displacement vector. 
The offset of 20 pixels in y-direction entered in the simulations 
is easily recovered as the location of the peak in the case of 
the 45-pix filter. For the maps filtered on a 5 pixel scale, the 
uncorrelated noise structures are strong so that one cannot un¬ 
ambiguously distinguish the maximum of the WWCC at the 
correct offset from other maxima caused by the noise. As a 
consequence, the position of the correlation peak can not be re¬ 
liably determined in this case. This prevents the determination 
of the displacement vector on noise-dominated scales less of 
than about 8 pixels. 

Cross-correlation spectrum 
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(b) 

Fig. 5. Cross-correlation functions for the pairs of wavelet- 
filtered maps from Fig.O The top panel (a) shows C w (t, 5pix), 
the bottom panel (b) C w (t, 45pix). The red crosses denote the 
center of the offset plane t - 0, where one can read the wavelet 
cross-correlation coefficient r(/) = C w (t = 0, /). The peak of 
the cross-correlation functions is displaced from the center by 
r(/) = Ti as the WWCC recovers the input displacement of the 
structure in the second map. 


The CC coefficients are extracted from the central pixel of 
the WWCC function (Eq. (fl9l) , crosses in Fig. 0) for all possible 
wavelet-filter sizes, /. In Fig. [6] (a,b) we show the dependence 
of the CC coefficient r(J) on the scale for the circular structures 
with S/N = oo and S/N = 5. As the two circles in the example 
are clearly adjacent, the cross correlation is negative at small 
scales, showing the expected anticorrelation below the mutual 
shift of 20 pixels and up to ~ 23 pix. The correlation curve has 
a shallow minimum at« 15 pix and monotonically increases up 
to 0.6 at large scales where both structures are strongly blurred 
through the wavelet filter, so that they turn statistically similar. 
As the offset is comparable to the diameter of the smaller circle 
/ p , the coefficient remains well below unity even at large scales. 
In case of S/N = oo, the anti-correlation is detected at all small 
scales, in the case of S/N = 5 (Fig. 0(b)) the correlation van¬ 



0 fi M IS 40' SO §§ 



0 10" 25 30 40 SO: ©5 


Fig. 6. Correlation coefficient r = C(r = 0) as a function 
of scale for the maps from Fig. [T] (a) represents the pure 
Gaussians, (b) the maps with noise resulting in S/N = 5. 

ishes at scales / & 5 pix due to the uncorrelated small-scale 
noise in the two maps. 

There are two main parameters shaping r(/), one is the dis¬ 
placement between the structures and one is the ratio of promi¬ 
nent scales of these structures, To show the effect of 

the displacement on the correlation coefficient r(/), we com¬ 
pute the WWCC for pairs of circular structures (cr\ = 10 pix 
and <j\ — 5 pix) mutually displaced by offsets of 0, 5, 10, 15, 
and 20 pix (Fig. [7]). 

The correlation obviously becomes stronger on all scales 
if the displacement between the structures is smaller. An anti¬ 
correlation occurs whenever if the mutual shift of the two struc¬ 
tures is larger than 1/3 of the size of the bigger one, almost 
independent of the size of the smaller Gaussian. For smaller 
displacements, the correlation remains positive at small scales. 
Without anticorrelation we always find a monotonic increase 
to large scales. When the offset between circular structures in¬ 
creases the CC coefficient is reduced on all scales and turns to 
negative on small scales. 

To quantify the behaviour of the correlation on small scales 
we examine the changes of the CC coefficient at the smallest 
wavelet-filter size, r(2 pix), as a function of the size ratio of the 
circles and their normalised displacement in Fig. [8j The struc¬ 
tures correlate strongly if their prominent sizes are comparable 
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Fig. 7. Correlation coefficient as a function of scale for two cir¬ 
cular structures (<xi = 10 pix and cr 2 = 5 pix) shifted by 0, 
5, 10, 15, and 20 pix (full line, dotted, dashed, dash-dotted, 
and dash-dot-dot-doted lines, respectively). The dashed vertical 
lines denote the prominent scale of the two circular structures. 



r/i 


Fig. 8. Contour diagram of the small-scale CC coefficient, 
r(2pix), as a function of prominent size ratio of circu] 
structures (/ P 2 // P i) and displacement normalised to the largest 
prominent scale (t;// p i). Here we used S/N = oo and l v \ = 60 
pix. 


Fig. 9. Contour diagram of the normalised cross correlation 
root /cco/^pi as a function of prominent size ratio of circular 
structures (/ P 2 // P i) and displacement normalised to the largest 
prominent scale (r;// p i). S/N = oo and / p i = 60 pix. 
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Fig. 10. Displacement vector as a function of scale for the pairs 
of Gaussian maps from Fig. \T\ (a) represents the results for 
pure Gaus^ians, (b) for the maps with S/N = 5. The vectors are 
shown in (he x - y-plane for individual discrete lags indicated 
on the x-a(ds at the origin of the plotted vectors. 


and the offset between them is small. The correlation weakens 
with decreasing the size ratio or increasing the offset. The cor¬ 
relation is about zero along the line / P 2 // P i ~ 0.48r;// p i - 0.21 
and it turns to negative for larger displacements reaching a min¬ 
imum at Ti/l p i ~ 0.5 and / P 2 // P i ~ 1. 

We can also characterize the anticorrelation due to dis¬ 
placement through the parameter Iqc 0 , defined as the scale at 
which the correlation coefficient r(l) turns from negative to 
positive, i.e. r(/cc 0 ) = 0- In Fig. 0 we show the dependence 
of /cc 0 on the relative displacement T;// p \ and the size ratio of 
the Gaussians / P 2 // P i. The scale of the cross-correlation root 
depends strongly on the offset, but hardly on the size ratio of 
the structures. For offsets below r,-// p i « 0.3 no root is found 
as the cross-correlation function is positive even at small fil¬ 


ter sizes. [For offsets larger than « 0.8/ p i, i.e. displacements 
comparable to the prominent scale of the largest structure, the 
cross-correlation function remains negative on all scales. At 
Ti * 0.7/ p i we find / C c 0 ~ l v i- 


The spectrum of displacement vectors 

Finally, we examine the spectrum of displacement vectors, 
r(/), measured for all wavelet-filter sizes. The displacement 
vector gives the coordinate of the peak of the WWCC rela¬ 
tive to the central pixel (Eq. (l20l)). Fig. [T0| shows the recov¬ 
ered displacement vectors between circular Gaussian structures 
for S/N = oo and S/N = 5. In case of infinite signal-to-noise 
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Fig. 11. Upper limit for the displacement scale vs. prominent 
scale estimated for seven pairs of equal size circular structures 
(having cr = 2.5, 5, 10, 15, 20, 25 and 30 pix) offset by r y = 4, 
12, and 20 pix. 

(Fig. |T0l(a)), the vectors are accurately recovered at scales less 
than the / p i ^40 pix. At larger scales the displacement vector 
is slightly underestimated by up to two pixels. This is due to 
the finite map size truncating parts of the filtered circles that 
fall beyond the map boundaries. More and more of the filtered 
structure shows up beyond the map edges when going to larger 
wavelet-filter sizes (see Fig. |2]) so that for the remaining part 
within the map boundaries the mutual shift appears smaller 
than the applied offset. The coinciding finite map boundaries of 
both maps impose a correlation for larger wavelet-filter sizes. 
This mimics a somewhat lower displacement of the structures 
and limits the applicability of the WWCC at large scales (see 
below) thereby affecting maps with low spatial dynamic range. 

The spectrum of the displacement vectors for noisy maps 
(S/N = 5; Fig.flQl(b)) shows completely random results for the 
noise dominated scales below 6 pixels where the correlation co¬ 
efficient is close to zero. Up to the minimum of the A-variance 
spectrum at 8 pixels, the displacement vectors remain unreli¬ 
able. At scales between 8 pix and 12 pix, where the noise is 
still strong, but not dominant any more, we can recover the dis¬ 
placement vector with an accuracy < 15 %. Up to 12 pixels, the 
cross correlation coefficient is also dominated by noise, rather 
than actual structures in the map (see Fig. 0(b)). For all scales 
above 12 pixels, the displacement vector is recovered with an 
accuracy better than 10 %. 

While the finite maps size poses an upper limit to the usable 
scales, the observational noise thus poses a lower limit to the 
scales for which we can determine the correct structure offset. 

4.1.4. Limiting scales 

To systematically quantify the limits for the accurate determi¬ 
nation of the spectrum of displacement vectors we study the 
minimum /i ow and maximum / upp scales for which the displace¬ 


ment vector can be recovered with an accuracy better than 

10 %. 

The upper limit for the scale, / upp 

Here, we consider the edge effects for different sizes of 
circular structures and different positions near the edge of the 
map. We use seven pairs of equal size Gaussian clouds simu¬ 
lated with S/N = oo and cr = 2.5, 5, 10, 15, 20, 25, and 30 pix 
(/ p = 8, 19,40, 60, 80, 100, and 120 pix). We displace each pair 
of maps in y direction, by T; = 4,12, and 20 pix thus moving the 
structures closer to the map edges. For every given offset r /9 we 
compute the spectra of displacement vectors for the seven pairs 
and determine the upper limit / upp where the measured spatial 
displacement r is recovered with an accuracy better than 10 %. 
This is performed by evaluating a circle of radius Sr = 0.1 x |t;| 
around the input displacement vector and testing whether the 
recovered vector falls into that circle. 

The maximum scale / upp at which the displacement vector 
can be recovered with 10 % accuracy depends on the truncation 
of the wavelet-filtered circles at the map edges. In Fig. [IT] we 
show the measured limiting scale / upp for a number of offsets 
and prominent scale sizes. We find a general decrease of the 
limiting scale / upp with increasing prominent structure size. For 
structures with prominent scales < 1 / 6 of the map size (/ p < 20 
pix) the displacement vector can be accurately measured up to 
scale / upp of about half the map size. For structures with larger 
prominent scales (> 50 pix) a reliable measurement of the dis¬ 
placement is only possible on smaller scales. The boundary 
effects strongly limit the usable scales. The truncation of the 
structure itself by the finite size of the map will result in an 
underestimation of the displacement at large scales. Circular 
structures with sizes comparable to the map size (> 100 pix) 
already fall beyond the map boundaries when offset by 12-20 
pix, i.e. the offset can only be measured in a small range of 
scales / upp ^ 10 pix. However, these are scales that may be 
dominated by noise already. 

To guarantee a reliable determination of the displacement, 
we fit the lower limit of the upper scale / upp as a function of the 
prominent scale / p by the linear relation 

4pp = -0.38 / p + 64pix = -0.38 l v + / Mp /2. (27) 

shown as dashed line in Fig. |TT] This relation can be used to 
constrain the / upp of circular structures whenever the displace¬ 
ment does not shift the prominent structures off the map bound¬ 
aries. 

The lower limit for the scale, /i ow 

In Fig.[l0](b) we have seen that for the S/N = 5 example the 
CC coefficient and displacement vector can not be recovered 
at scales below the scale of the minimum of the A-variance 
spectrum / m i n = 8 pix (Fig. iu>. 

Figure [12] shows the lowest scale /i ow for which the dis- 
placement vector can be retrieved with an accuracy better than 
10% as a function of the scale / m i n for a number of different 
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Fig. 12. Lower limit for the displacement scale (/i ow ) as a 
function of the A-variance minimum scale (/ m in) for S/N = 
2.5,5,10,25, and 50, and displacements of r t = 4,12,20 pix 
(black, gray, and open symbols). An upper limit to the lower 
scale is described by the solid line: how = 1.5 / m in + 1.5. 


Fig. 13. Measured mean displacement ((r)) vs. given displace¬ 
ment (rO for S/N = oo (circles), S/N = 25 (squares), and 
S/N = 5 (triangles). The circles and triangles are displaced by 
1 pixel from the identity line (dashed line) for a better visibility. 
The error bars represent 1 cr error of the mean offset. 


noise levels and displacements, using equal size circular struc¬ 
tures. We find a considerable scatter of the limit for the reliable 
determination of the displacement vector when considering dif¬ 
ferent offsets, but we can give a safe upper limit to how by 

ho w = 1.5 (/ min + !)• (28) 


offsets, i.e. we can reliably determine the displacement vector 
between two structures t(Z) in the full range between l\ow and 

Ajpp- 

4.2. Structures of fractional Brownian motion 


shown as solid line in Fig.fl2l 

The range l m m ^ l ^ how thus characterizes scales where 
the noise does no longer hide the underlying structure - as for 
scales l < /min - but where it still significantly influences the 
cross-correlation spectra and the measured displacement vec¬ 
tors. 

To finally verify the reliability of the recovery of the dis¬ 
placement vector within the limits given by Eqs. (l27l) and (l28k 
we restrict the spectrum of displacement vectors to that range 
when evaluating the results for 315 pairs of maps covering 
all combinations of Gaussian widths cr = 1, 4, 7, 10, 13, and 
16 pix (Z p = 4, 16, 28, 40, 52, and 64 pix), input offsets of 
Ti t y = 0,5,10,15, and 20 pix, and S/N values of oo, 25, and 
5. Fig. [T3l shows the measured mean displacement (r) against 
the given displacement r* for the ensembles of 21 pairs hav¬ 
ing the same offset and noise level. The standard deviation of 
the derived displacements is very small so that the plotted error 
bars basically fall within the size of the symbols so that they 
are hard to be seen. 

Overall we find a very good agreement with a minimal scat¬ 
ter, proving the method to work, but a slightly overestimated 
displacement at r* = 0 pix for high noise (S/N = 5, triangle) 
and a slightly underestimated displacement for the largest off¬ 
sets at all noise levels, due to the finite map size effects. The 
deviations fall, however, well below one pixel size. Hence, we 
conclude that the cutoff of the spectrum of the displacement 
vectors by low and upper limits allows us to recover accurate 


In contrast to the circular Gaussian structures that have a single 
well-defined prominent scale, maps of the interstellar medium 
often show a self-similar behavior without individual promi¬ 
nent scales. They can be represented by the fractal structure of 
fractional Brownian motion maps. They are characterized by a 
power-law distribution of scales and are periodic by generation 
in Fourier space. We use a spectral index f = 3 representing 
a typical me an v alue observed in the cold interstellar medium 
( Falgarone et al.ll2QQ7 ). As for the circular Gaussian structures 
we evaluate how well the WWCC analysis can detect scales en¬ 
hanced in one structure compared to the other and the displace¬ 
ment of structures at particular scales within the two maps. 


4.2.1. fBm structures displaced on different scales 

In real clouds, structures can be offset with respect to each 
other on different spatial scales. This occurs, e.g., in differ¬ 
ent velocity channel maps of line observations where the large 
scale structure is affected by systematic motions but the small 
scale structure by turbulent motions, in dust emission maps at 
different wavelengths where local temperature gradients cre¬ 
ate small scale displacements independent of the possible large 
scale structure, or in maps of different chemical tracers where 
incident UV radiation from one direction creates a chemical 
gradient on scales not mixed by turbulent flows. 

To test this, we start from an fBm map with a spectral index 
£ = 3 and S/N = oo (Fig.[l4l left panel), and shift all structures 
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Fig. 14. fBm map with £ = 3 (left panel) and the same map 
after displacing all structures on scales l > 10 pix by = 10 
pix and r ity = -15 pix (right panel). This pair is to mimic an 
interstellar cloud structure that could be created by the super¬ 
position of small scale turbulence and a large scale velocity 
gradient. 



Fig. 15. A-variance spectra for the maps from Fig.fl4l 

with scales larger than 10 pixels by |r z -| = 18 pix (r i;JC = 10 pix 
and Tiy = -15 pix) using the Fourier shift theorem (Eq. (l25l)). 
The resulting map is shown in Fig. |T4| (right panel). The spatial 
shift of structures does not affect the relative contribution of 
structures as a function of their size, so that both maps have an 
identical A-variance spectrum, shown in Fig. Q3] It is charac¬ 
terized by a perfect power-law with an exponent £-2= 1, rep¬ 
resenting the fully self-similar scaling up to the largest scales 
where the structures are limited by the available map size. 

Figure [16| shows the results of the WWCC for the com¬ 
parison of both maps. The correlation coefficient is unity at 
small scales (because the structures are not offset) and gradu¬ 
ally decreases to a minimum at 17-18 pix due to displacement 
by 18 pix between the structures (Fig.[l6|(a)). On large scales 
the structures are blurred through the wavelet filter to sizes that 
exceed the displacement amplitude. This leads to monotonic 
increase of the correlation there. The displacement vector is 
perfectly recovered for all scales, being zero on scales / < 10 
pix and r x = 10 pix, r y = -15 pix for larger scales (Fig.[T6l(b)). 

Fig. U3 shows the corresponding results when adding noise 
to the maps. The noise level of S/N = 5 appears as an increased 
A-variance at small scales (/ < 8 pix, panel (b)). It lowers the 
CC spectrum on scales l < 8 pix, but has no impact on the mea¬ 
surement of the displacement vector. As the map is filled with 
structures on all scales, the impact of the noise is much lower 



0 10 2Q 30. 40 SO 60 



Fig. 16. Results of the WWCC for the pair of maps from 
Fig. |T4j (a) correlation coefficient and (d) displacement vector 
as a function of scale. 1 cr error bars are given. 

here than for the sparsely populated maps inhabiting a single 
Gaussian circular structure from Sect. 14.11 To ensure that we 
actually find a smooth transition between the two extreme cases 
of the single Gaussian circular structure and the fBm map, we 
also checked the impact of the noise on recovery of the dis¬ 
placement spectrum in case of multiple Gaussian structures. 
When we compare two maps each having three randomly lo¬ 
cated Gaussian circular structures with the same characteristics 
as in Fig.|T|(b) and measure the lower limit for the displacement 
scale we obtain /i ow = 5 pix (/ m i n = 7.7 pix), i.e. a value that 
is considerably smaller than that for the single Gaussian struc¬ 
tures (/i ow = 12 pix, see open upward triangle in Fig. [12]). For 
pairs of maps inhabiting 30 Gaussian structures the impact of 
the noise further decreases. We obtain /i ow = 4 pix (/ m i n = 6.3 
pix) confirming that the noise dominates smaller and smaller 
scales if the maps are filled with more circular structures. This 
is in line with a negligible noise impact for the fBm maps filled 
with structures on all scales (Figs.[l6l(b) andfT71(d)). 

Identification of correlated structures 

To understand the nature of the correlations quantified by 
the WWCC, it is useful to visually inspect the filtered maps 
for the individual filter scales /. In the original maps (Fig. [17] 
(a)) it is almost impossible to locate individual structures that 
are correlated between the maps. In the corresponding maps 
that were wavelet-filtered on the scale of 6 pix (Fig. [T8ltop). 
the scale where the CC spectrum showed a peak indicating 
many correlated structures (r(/ = 6 pix) « 0.6, Fig. [17] (c)), it 
may be possible to visually identify some prominent structures, 
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Fig. 18. Top left. The fBm map (Fig. [T7](a), left) filtered on 
scale of 6 pix. Top right. The fBm map displaced on scales 
> 10 pix (Fig. [171(a). right) filtered on scale of 6 pix. Bottom 
left. The product map generated by multiplying the intensities 
of the top left and top right fBm maps at each pixel (Eq. (28)). 



(c) 


L (pixel) 



Fig. 19. Top left. The fBm map (Fig. [T7](a), left) filtered on 
scale of 17 pix. Top right. The fBm map displaced on scales 
>10 pix by Ti x = -10 pix and r^ y = 15 pix (Fig. [171(a). right) 
filtered on scale of 17 pix. Bottom left. The product map of two 
m maps (top left and top right). Bottom right. The product 
p of the top left and top right maps where the latter is re¬ 
versely shifted by r x = -10 pix and r y = 15 pix (Eq. (29)). 


(d) 
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Fig. 17. (a) Original fBm map generated with S/N = 5 (left) 
and the same fBm map displaced by r itX =10 pix and r^ y = 
-15 pix on scales l > 10 pix (right), (b) A-variance of the orig¬ 
inal and displaced fBm maps, (c) correlation coefficient and 
(d) displacement vector as a function of scale, lcr error bars 
are presented. 


but this is extremely difficult due to the noise contribution. An 
easy identification of correlated structures is however possible 
if we visualize the integration kernel of the correlation function 
Eq.® for zero offsets r = 0, i.e. show the product of the two 


maps at each pixel, 

P(x , /) = F(x , /) x G(x, /). (29) 

In the product map P(x , /) (Fig. UH1 bottom panel) only rel¬ 
atively strong (positive or negative) features that agree between 
both maps show up@. One can easily recognize individual small 
“clumps” present in both maps at the same location that can be 
identified by comparing the product map with the original maps 

(Fig.[I7](a)). 

In Fig. [19] we repeat the experiment for the larger scale of 
/ = 17 pix where the correlation coefficient is at minimum 

2 In observed intensity maps, there are of course no significant neg¬ 
ative structures, so that we should only see the bright structures there. 
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(Fig. H7](c)) because the structures are displaced by r^ x =10 
pix and r^ y = -15 pix (|t; = 18 pix, Fig.[T7](d)). The product 
map of the filtered maps (Eq. d29lh Fig. [I9l bottom left panel) 
shows no structures that can be identified in the individual fil¬ 
tered maps or the original maps in Fig.[T7](a) and the amplitude 
of the structures in the product map is relatively small. When 
using, however, the recovered displacement vector at l = 17 pix 
to shift the second map back by the measured offset vector r(/) 
and take the product of the shifted map G(x + r(/), /) and the 
first map F(x , /) (top right), 

P{x , t(Z), 0 = F(x, l) x G(x + t(Z), /), (30) 

we find that the product map (Fig. [191 bottom right) recovers 
all of the relatively strong (positive or negative) structures seen 
in the filtered maps. This allows us to localize individual corre¬ 
lated structures to address their shape and origin, independent 
of a mutual shift. In Sect. [5] we encounter a pair of observed 
maps of the G 333 molecular cloud with properties that are very 
similar to the example considered here. In the observed maps 
we also find correlated structures with matching locations at 
small scales, but a global displacements when considering large 
scale structures. 

4.3. Structures of fBm clouds with enhanced scales 

Pronounced scales in otherwise self-similar clouds may indi¬ 
cate special physical processes acting on those scales, therefore 
giving access to important characteristics of observed maps. 
Being able to find and compare those scales in different maps 
is therefore an essential step in understanding interstellar tur¬ 
bulence. One process producing such enhanced scales is the 
opacity of partially optically thick lines. They provide a satu¬ 
rated picture for column densities above a specific threshold. 
To mimic this opacity saturation effect in simulations, we en- 



Fig. 20. Top left. A pure fBm map generated with S/N = oo and 
spectral index £ = 3. Top right. A pure fBm map filtered with 
maximum filter of size 15 pix. Bottom left. efBm map which is 
generated by superimposing the pure fBm structure (top left) 
and the filtered fBm structure (top right) and for brightness pa¬ 
rameter W = 1. Bottom right. efBm map shifted by r^ x =10 
pix and r^ y = 10 pix. 



Fig. 21. A-variance spectra of the pure fBm structure and fBm 
structure filtered at the scale of 15 pix by the maximum filter 
(full and dotted lines, respectively). 

hance spatial structures in an fBm map F(x) by smoothing it 
with a maximum filter described by 

MifF(x )) = max x' : \x' - x\ < lf/2F(x') , (31) 

where If is the filter size. The maximum filter enhances the 
signal on the scale of the filter size and washes out structure 
on scales less than the filter size. In this way we introduce a 
prominent scale with size If into the map. This is demonstrated 
in Fig. [20l The upper left panel shows the original fBm struc¬ 
ture and the upper right panel the result of the convolution with 
a maximum filter of size If = 15 pix. Structures below If pix 
are washed there while circles of 15 pix diameter represent the 
dominant structures now. As the pure opacity effect is always 
superimposed by additional small scale variations this would 
give an unrealistic picture. We construct a more realistic picture 
by combining the original fBm, F(x), with the scale-enhanced 
fBm map, Mi f (F(x))\ 

F m (x) = WF{x) + (2 - W)M k {F(x)\ (32) 

where we introduce the inverse brightness contrast parame¬ 
ter W, (0 < W < 2), characterizing the contribution of the 
pure fBm in the total map. Then, 2 - W is the relative bright¬ 
ness of the scale-enhanced map. For values of W < 1 the en¬ 
hanced structure is brighter than the pure fBm structure while 
for W > 1 it is fainter. This combination should account for the 
characteristics of the globally self-similar structure and opacity 
effects. The resulting fBm with enhanced scale (in the follow¬ 
ing efBm) map is shown in the lower left panel of Fig.l20lfor 
W= 1. 

4.3.1. Recovering the enhanced scale. 

In our numerical experiment we know the size of the maximum 
filter, but in real observations, one first has to detect and quan¬ 
tify the prominent structures, i.e. the manipulation of the data in 
our case. [Mac Low & Ossenkop f ( 2000 ) and 008 have shown 
that the A-variance is an appropriate tool to detect these scales 
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Fig. 22. A-variance spectrum cr 2 A (l) of a filled circle of diameter 
20 pix (full line) and derivative of the logarithmic A-variance 
spectrum g(l) of the filled circle (dotted line) for S/N = oo. 
Vertical lines denote the measured prominent scale of the filled 
circle at / p = 21.9 pix and the critical scale l c = 10.3 pix char¬ 
acterizing the peak of the gradient spectrum (Eq. (I34l0. 


in general astrophysical data sets with irregular boundaries and 
variable data reliability across the maps. The A-variance spec¬ 
trum measures all prominent scales that are strong enough to 
create a peak in the spectrum. In a globally self-similar struc¬ 
ture like our fBm (or efBm) maps, there is, however, a mono- 
tonically increasing A-variance spectrum. Figure I2T1 compares 
the A-variance spectra of the original fBm with the spectrum 
obtained for the efBm map from Fig. [20] We find no peak in the 
A-variance spectrum - the underlying fBm structure dominates 
- but information about the filter size, If , is apparently present 
in the relative drop of structure variation at scales below the 
filter size. The convolution of the map with the maximum filter 
of size If = 15 pix filters out smaller scale structure. The efBm 
has less structure on scales <15 pix, while at scales larger than 
the size of the maximum filter the A-variance of the original 
fBm is basically recovered in the efBm. 

Instead of searching for a peak in the A-variance spectrum, 
we therefore need a new approach to find enhanced scales by 
looking at the steepening of the A-variance spectrum that is 
visible in Fig.|2T] As an auxiliary step, we first look at the pure 
maximum filter that is described by a filled circle with given 
diameter If. 

Circular structures with constant intensity 

The A-variance spectrum of a filled circle of diameter d = 
20 pix is shown in Fig. [22] as a full line. The maximum of the 
A-variance falls at / p ~ 20 pix « / f , as predicted by 008. We 
verified this for filled circles with If = 4, 8,12,16,20,24,28 
pix. 

To measure the prominent scale even without clear A- 
variance peak, we examine the gradient of the logarithm of the 


Fig. 23. Relation between critical scales and diameter of the 
circles for S/N = oo, 25, and 5. The best-fit linear regression 
line (/ c « 0.52 / p ) for S/N = oo is shown as solid line. 


A-variance spectrum, 


g(Q = 


d log o-2(/) 

dlog/ 


(33) 


in Fig. [22] (dotted line). The steepening of the A-variance spec¬ 
trum below the prominent scale If is reflected by a pronounced 
peak in the gradient spectrum. The root g = 0 corresponds to 
the prominent scale / p . We define the critical scale as the max¬ 
imum of the gradient spectrum: 


l c = arg max 


dlog 0-2 (Z) 
dlogZ 


(34) 


For the filled circle with 20 pixels diameter the measured crit¬ 
ical scale is l c = 10.3 pix, approximately half of the circle’s 
diameter. 

When computing the A-variance gradient spectra for circles 
of If = 4, 8,12,16,20,24,28 pix (circles in Fig.l23lk we find a 
linear relation between measured l c and If 


l c = (0.52 ± 0.02) Z p - (0.14 ± 0.36) * 0.52 Z f , (35) 
shown as solid line in Fig. [23] 


Noisy data 


To test the robustness of this approach against observational 
noise, we add different levels of white noise to the circle maps 
and compute their A-variance gradient spectra. In noisy maps, 
the A-variance spectrum at small scales is dominated by noise 
fluctuations (compare Figs. [4] and [3]) so that the critical scale 
l c is not directly measurable through Eq. (l34l) . To still derive 
Z c , one needs to subtract the noise contribution from the A- 
variance spectra of the contaminated circular structures: 


^a(0si 


— a (Omap \(J)n 


(36) 


The corresponding noise spectrum can be obtained by running 
the A-variance analysis on the noise map, obtained e.g. from 
emission-free channels in an observed spectral line cube. As 
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Fig. 24. Mean critical scale l c measured for eight fBm maps 
filtered with the maximum filter of size If for S/N = oo, 25 and 
5 (circles, squares, and triangles, respectively). 1 cr error bars 
measured in a set of 10 different realizations are superimposed. 


we construct the noisy map in our numerical simulations by 
simply adding a noise map to the original map we do not need 
to extract it separately here. 

Using this noise subtraction, we repeated the computa¬ 
tion of the gradient peak scale / c from maps with different 
circle sizes and noise levels S/N = 25 (squares in Fig. l23l) 
and S/N = 5 (triangles). We find a very good match to the 
noise-free results with deviations by up to at most 2 pixels 
induced by the noise in the data. The regression coefficient 
(«a = 0.54 ± 0.02) for the noisy maps coincides with one for 
S/N = oo (a = 0.52 ± 0.02) within the error limits. 

For observed maps the determination of the noise A- 
variance <x^(/) no ise may, however, be affected by uncertainties 
that will propagate into the determination of the critical scale 
l c if the latter falls into the noise-dominated regime. Errors in 
the noise subtraction will affect the A-variance gradient spec¬ 
trum at scales below the A-variance minimum, / < l m i n (see e.g. 
Fig- 021(b), /noisedominated £ /min = 5 pix). This puts a lower limit 
on the scale l > l c = / min that can be reliably measured. 

efBm structures 

In Figure [24] we show the result of the equivalent experi¬ 
ments for the fBm structures with enhanced scales and an in¬ 
verse brightness contrast W = 1. The efBms are generated from 
a fBm map with spectral index £ = 3 which is convolved with 
the maximum filter of the sizes If = 5, 9, 13, 17, 21, 25, 29 pix. 

For each efBm we use the gradient of the A-variance spec¬ 
trum to measure the critical scale (Eq. (l34l)). Measured critical 
scales against maximum-filter sizes are plotted in Fig. [24] The 
mean critical scale and its uncertainty are estimated from 10 
different random realizations of efBms. The circles with error 
bars represent the noise-free results. The data can be described 


by a linear relation, 


<Z C > = (0.53 ± 0.01) k - (0.61 ± 0.1), 


(37) 


represented by the dashed line in Fig. [24] Within 1 pix accuracy 
this is identical to the simple approximation without offset from 
the origin l c « 0.52 If (in agreement with Eq. (135th obtained for 
the circular filter, shown as solid line in Fig. [24] 

The squares and triangles in Fig. [24] show the measured 
critical scales for the equivalent noisy maps with S/N = 25 
and S/N = 5 after applying the noise correction from Eq. (l36t . 
One can see that the squares and triangles coincide well with 
circles for If > 1 pix and they can be fit by the same linear 
relation given by Eq. <03 and <03 (see Fig. El. 

To verify that the relation is robust against changes of the 
efBm spectral index and the inverse brightness contrast W we 
repeated the experiment for a range of parameters. We vary the 
fBm spectral index in the rang e of 2.5 to 3.5 typical for inter¬ 
stellar clouds dFalgarone et al.l2007l) and brightness parameters 
W = 0.1,0.25,0.5,1,1.5,1.75,1.9. The previous experiment 
used £ = 3 and W = 1. For each combination of £ and W , we 
recover the critical scales l c for If = 5, 9, 13, 17, 21, 25, 29 
pix and compute the slope a of the / c - If relation in the range 
If > / m in pix. To estimate the statistical uncertainty of a we used 
10 different random realizations so that the error is composed 
from the ensemble variation and the imperfection of the linear 
fit. Figure [25] shows the results for S/N = oo and S/N = 5. In 
case of S/N = oo (upper panel), the slope varies in the range 
between 0.47 and 0.63. It tends to increase with decreasing the 
spectral index £ and increasing inverse brightness contrast W. 
Averaged over the full range of spectral indexes (£ = 2.5-3.5) 
and inverse brightness contrast parameters (W = 0.1 - 1.9) we 
obtain a mean slope of (a) = 0.54 + 0.03. When knowing £ 
from the overall slope of the A-variance spectrum, the gradient 
can be even better constrained. For the three spectral indexes 
we obtain (a)^= 2.5 = 0.57 ± 0.03, (a )^=3 = 0.54 ± 0.02, and 
(a)£= 3 £ = 0.51 ± 0.01. These values can be used directly to 
measure the enhanced scale from the critical scale in the A- 
variance spectrum for large S/N ratios. 

In case of the noisy maps (lower panel in Fig. [25]), the vari¬ 
ation of the slope is only slightly larger for inverse brightness 
contrasts W < 1.5 but strongly deviating from the noise-free 
behaviour for low contrasts of the enhanced structure, i.e. val¬ 
ues of W > 1.75. This is intuitively clear as we see less and 
less enhanced structure in a map if we reduce its contribution 
by increasing W to values close to two. If the enhanced struc¬ 
ture is strongly diluted (W > 1.5), the efBm maps appear more 
like pure fBm maps with only small variations of the gradient 
of the log cr^(/) so that the determination of the gradient peak 
/ c , becomes very sensitive to the noise contribution. Additional 
noise makes a reliable detection of the enhanced scale increas¬ 
ingly difficult. 

The impact of the inverse brightness contrast on the A- 
variance gradient is visualized for the noise-free case (S/N = 
00 ) in Fig. [26] It shows the gradient spectra of four efBm maps 
with £ = 3, Zf = 15 pix, and inverse brightness contrasts 
W = 0.25,1,1.75,2. The case of W = 2 (dot-dashed line) 
represents the pure fBm structure showing the co nstant gra di- 
ent g(l) = £ - 2 = 1 for two-dimensional maps (IStutzki et al 
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Fig. 25. Measured slope a against inverse brightness contrast 
W for efBm maps having spectral indexes £ = 2.5, 3, 3.5 (full, 
dotted, and dashed lines, respectively). The upper plot shows 
the result for S/N = oo, the lower plot for S/N = 5. The data 
points for f = 3 and 3.5 (squares and triangles) are horizontally 
displaced by W = 0.02 and W = 0.04, respectively, from the 
£ = 2.5 data points (diamonds) for a better visibility. The error 
bars represent 1 cr variations in the statistical sample. 

19981) at scales up to 30 pixels (dashed vertical line in Fig. [26]). 
The drop towards larger scales is due to the impact of the map 
boundary limiting to size of large structures (008). This is con¬ 
sistent with our results from Sect. 14.1.41 that only prominent 
scales up to 64 pixels, i.e. half the map size can be detected. 
To omit the boundary effects we restrict the analysis to criti¬ 
cal scales l < 30 pix. When adding scale-enhanced structures, 
W < 2, we see a growing peak in the gradient spectrum at 
l c = 8 pix. It turns sharper if we lower W, i.e. the contrast of 
the peak gradient g max relative to the average gradient (g) grows 
when we increase the contribution of the scale-enhanced struc¬ 
tures by lowering the inverse brightness contrast parameter W. 
The contrast of the peak value g max relative to the average A- 
variance gradient (g) determines whether l c can be reliably de¬ 
termined, even if uncertainties are added to the spectrum in 
terms of noise. 

In Fig. [27] we show the result of the systematic study of 
the relation between g max and (g) for different values of the 



Fig. 26. Gradient spectra of fBm maps generated with £ = 3, 
S/N = oo, enhanced on scale of If = 15 pix and for differ¬ 
ent inverse brightness contrast parameters W = 0.25,1,1.75 
and 2 (full, dotted, dashed, and dot-dashed lines, respectively). 
The vertical dashed line indicates the upper limit to the scale 
where the A-variance slope can be measured reliably. Above 
that scale, the spectrum is strongly affected by boundary effects 
of the assumed finite map. 



Fig. 27. S/N = 5. Measured maximum gradient versus mean 
gradient for efBm maps having spectral indexes £ = 2.5, 3, 3.5 
(full line, dotted and dashed lines, respectively) and different 
inverse brightness contrasts W = 0.1,0.25,0.5,1,1.5,1.75,1.9. 
The error bars indicate the standard deviation of the ensem¬ 
ble of 10 different random realizations and seven different en¬ 
hancement scales. The long-dashed line is the identity repre¬ 
senting straight spectra, i.e. g max = (g). The dash-dot line 
represents the lower limit of g max as a function of (g) below 
which the scale of the peak cannot be reliably determined (see 
Eq. ©). 


spectral index £ and the inverse brightness contrast W for 
S/N = 5 where (g) is measured between / m i n and 30 pixels. 
Averages and error bars are computed for efBms ensembles 
from ten different random numbers and enhanced on scales 
of If = 5,9,13,17,21,25,29 pix. The series of different in- 



































T.G. Arshakian and V. Ossenkopf: WWCC analysis of structural scaling in turbulent clouds 


17 


verse brightness contrasts W for the same spectral index £ are 
connected by solid (f = 2.5), dotted (£ = 3.0), and dashed 
(£ = 3.5) lines. The long-dashed line indicates the identity, i.e. 
the behaviour of pure fBms (W = 2) not showing any maxi¬ 
mum in the gradient (dot-dashed line in Fig. [26]). Weakly en¬ 
hanced fBm maps (W = 1.9) are positioned slightly above the 
identity. The deviation grows for stronger enhanced fBms, i.e. 
towards smaller inverse brightness contrasts W < 2. 

As Fig. [25] has shown that for noisy data, the scale of the 
steepest gradient cannot be reliably determined if the efBm is 
only weakly enhanced, so that the gradient peak is too shallow, 
we can translate Fig. [27] into an easily usable criterion for the 
contrast of g max relative to (g) that still allows for a reliable 
measurement of a and thus k through / c . For S/N = 5, points 
with W > 1.75 had to be excluded. This corresponds to the 
area below the dash-dotted line in Fig. [27] It can be described 
by gmax values that fall below 

gmax, limit = (1.02 ± 0.01)<g> + 0.10 + 3c^ max> = 1.02(g) + 0.15 

(38) 

This relation is computed by connecting the data points with 
W = 1.75 and adding the 3 cr error margin of <x<g max ) = 0.017 
of the individual data points. All efBm structures where the 
enhancement scale cannot be measured reliably due to a too 
high inverse brightness contrast (W > 1.75 in Fig. [23(b)) fall 
below this relation. For W < 1.5 critical scales can be reliably 
recovered for S/N = 5. The mean slopes obtained for S/N = oo 
and S/N = 5 (0.54+0.03 and 0.51+0.04) agree within the error 
limits. 

Figure [27] can also be used to estimate spectral index and 
inverse brightness contrast of an efBm by measuring g max and 
(g). For high inverse brightness contrasts W > 0.5 correspond¬ 
ing to about g max - (g) < 0.4 spectral index and inverse bright¬ 
ness contrast can be uniquely determined, for lower values of 
W there is an ambiguity in determination of £ and W but it is 
still possible to determine f with an accuracy of about 0.5. 

One can generalise the approach of measuring the enhance¬ 
ment scale If in efBm maps (Eq. (l38l ) ) for S/N levels larger than 
S/N = 5. They show different thresholds of the inverse bright¬ 
ness contrast above which If can not be determined reliably 
anymore. When lowering the noise, the threshold approaches 
W = 2 approximately linearly with the noise level, shifting 
down the dot-dashed line in Fig. [26] towards the identity line. 
This corresponds to a change of the constant in Eq. (l38l ) from 
0.1 to 0 when S/N increases: 

gmax,limit - (g) + 0.5(S/N) 1 + 3 <X(g max ). (39) 

Global fit 

We can adopt the single value (a) = 0.52 + 0.04 for esti¬ 
mating the enhancement scale, 

k = lj( 0.52 + 0.04) = (1.92 + 0.15) Z c . (40) 

Even without a priori knowledge of the spectral index of an 
efBm cloud, one can measure its enhancement scale with accu¬ 
racy of about +8 % if the prominent structure is strong enough 
to meet the criterion of Eq. ([39b . 



Fig. 28. Spectrum of the CC coefficient for the original fBm 
map and efBm map from Fig. l20l (bottom left). 

However, Fig. [27] can also be used to estimate spectral in¬ 
dex, allowing for a more accurate determination of the slope 
a. When including noise levels from S/N = oo to S/N = 5 
in the fit, we measure the mean slopes (a)^= 2.5 = 0.55 + 0.04, 
(a)£=3 = 0.53 + 0.05, and (a)£= 3.5 = 0.49 + 0.06 allowing for a 
more accurate determination of the enhancement scale If than 
through the general equation (l40l) . 

Altogether this gives a recipe for the measurement of a pro¬ 
nounced scale in a structure that is globally dominated by self¬ 
similarity. After applying the noise correction to the A variance 
spectrum, the A-variance gradient spectrum has to be evalu¬ 
ated to measure the mean gradient (g), the peak gradient gmax, 
and the location of the peak / c . Depending on the noise in the 
data, the user has to evaluate whether the brightness contrast 
of the enhanced scale of a cloud is strong enough. If the gra¬ 
dient peak falls above Eq. ([39b . i.e. g max > gmax,limit, then it is 
sharp enough to allow for a reliable measurement of its posi¬ 
tion and one can use the critical scale l c > / m i n to recover the 
enhancement scale using the relation If = 1.92 l c (Eq. (|4Q b). 

4.3.2. Detecting and recovering displacements. 

Next, we explore the power of the WWCC to quantify the dis¬ 
placement between different efBm structures through the CC 
coefficient spectrum and the spectrum of displacement vectors. 

In Fig. [28] we show the CC coefficient spectrum when com¬ 
paring the original fBm and the efBm map without displace¬ 
ment (Fig. [20] top left and bottom left panels). The structures 
are highly correlated on small and large scales. The scale filter¬ 
ing with If = 15 pix produces a small minimum on a scale of 
about 5 pix due to the reduction of structures below the filter 
size. Above the filter size both maps turn very similar leading 
to a cross correlation coefficient close to unity. 

Figure[29](a) shows the result of the WWCC when compar¬ 
ing the fBm map and the efBm map displaced by r, = 14.1 pix 
(r ijX = 10 pix and r ity =10 pix, Fig. [20] top left and bottom 
right panels). We find a behaviour very similar to that for the 
displaced Gaussians in Fig. [6] The spectrum of the CC coef¬ 
ficients shows an anti-correlation at small scales with a mini- 
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Fig. 29. Results of the WWCC for the original fBm map and the 
shifted efBm map from Fig. [20] (bottom right): (a) correlation 
coefficient and (d) displacement vector as a function of scale, 
lcr error bars are given. 


mum at ~ 8 pix, a steep increase from the minimum to about 
30 pix, and a shallower increase at large scales. The spectrum 
turns to positive CC coefficients at about 14 pix approximately 
matching the offset amplitude. At larger scales, the correlation 
becomes stronger as the wavelet filter only measures structures 
with sizes above the enhancement and offset scale. It is evident 
that the displacement weakens the correlation on all scales as 
compared with Fig. [28] The effect is strongest at small scales, 
i.e. below the offset amplitude. The enhancement of the 15 pix 
scale in the efBm relative to the original fBm does not affect 
the recovery of the displacement vector (Fig.[29](b)). It is accu¬ 
rately measured over all scales l < 40 pix. At larger scales, the 
finite map size (see Sect. 14.1.31) leads to a small underestimate 
of the amplitude of the displacement vector by ~ 5 % < 1 pix. 

To combine all of this, we finally consider a pair of efBm 
maps with additional noise and an offset by t; = 14.1 pix 
( 'r iiX = 10 pix and = 10 pix, Fig. [30]). The first efBm 
map is generated using a maximum filter of size In = 15 pix, 
the second efBm map with maximum filter of size In = 5 
pix. The noise in both maps gives a S/N = 5. The maps are 
shown in Fig.[30](a), the A-variance spectra in[30](b), the cross¬ 
correlation spectrum in[30](c), and the spectrum of detected dis¬ 
placement vectors in[30](d). The effect of the noise is evident 
in the A-variance spectra as an increase of the variance towards 
small scales. Below / m i n ~ 6 pix the noise dominates, leading to 
a zero cross-correlation, but it has little effect on the recovery 
of the displacement vectors. The change of the vector ampli¬ 
tude remains within a 10 % accuracy. Above the noise scale 
the correlation curve behaves in the same way as in the case of 
S/N = oo (see Fig. [291(a)). The noise has no effect on r(l) at 


scales larger than « 6 pix. The amplitude of the displacement 
vector is also slightly underestimated on large scales as in the 
noise-free case. The scale enhancement in both structures has 
basically no effect on the outcome of the WWCC. It is only 
visible in the A-variance where the two critical scales can be 
measured as indicated in Fig. Eokbi. 

To study the quantitative impact of the displacement of the 
efBm maps on the CC coefficient more quantitatively, we per¬ 
form the same experiment as for the Gaussians in Fig. [7] using 
two efBm maps filtered with maximum-filter sizes of In - 15 
pix and In = 5 pix displacing them by r^ y = 0,5,10,15, and 
20 pix relative to each other. The result is shown in Fig. [31] The 
general behaviour is qualitatively similar to Fig. [7J The shape 
of the CC coefficient spectrum is regulated by the offset value 
and ratio between the smaller to the larger filter size, ln/k\ • 
Without any offset between the maps, we find a strong corre- 
tion between structures on all scales. It increases quickly to 
unity above the scale of the smallest filter size, In = 5. If the 
structures are offset by 5 pix (dotted line) there is no overlap 
any more below that scale so that the cross correlation coeffi¬ 
cient vanishes. The global shift leads to a lower correlation on 
all scales, including the largest ones. With increasing the offset 
(dashed, dashed-dotted, and dash-dot-dot-doted lines) the r(/) 
curve shifts to larger scales, a scale range of negative values 
(anticorrelation of filtered structures) appears, and the correla¬ 
tion becomes weaker on all scales. Obviously, the WWCC de¬ 
tects the displacement in the same way for the simple Gaussian 
structures in Fig. [7] and for close-to self-similar (e)fBm struc¬ 
tures. 

Finally, we test how well the displacement can be recov¬ 
ered for pairs of maps maximum-filtered on different scales. 
We filter the original fBm map with maximum filters of size 
k = 5,9,13,17,21,25,29 pix to generate seven efBm maps 
with £ = 3.0 and W = 1. Combining these seven maps (with¬ 
out repetition) gives 21 pairs of distinct maps. For a given pair 
of maps, one map is shifted by T; = 0,5,10,15,20 pix, respec¬ 
tively and for each shift we perform the WWCC and recover 
the displacement (r; Eq. (12Tb ) for ten random realizations of 
the fBm. Figure [32] shows the mean recovered displacements 
separately for S/N = oo, 25 and 5. 

Similar to the case of the Gaussian circular structures 
(Fig.ua, the mean displacements (r) are recovered with high 
accuracy for all noise levels and all offsets less than 1/6 of the 
map size (r, < 20 pix) also for the efBm structures. We also 
tested the identity between input and recovered displacement 
for efBm maps generated for all combinations of the spectral 
index ( = 2.5,3, and 3.5 and the inverse brightness contrast pa¬ 
rameter W = 0.5,1,1.5, (Eq. ([32l) ) and found a perfect agree¬ 
ment for all combinations of £ and W. The errors of the mea¬ 
sured offset depends on the parameter W. For large W the offset 
error slightly increases, however, it remains within the resolu¬ 
tion limit ±1 pix. 


3 In the ideal case considered here, a recovery of the critical scale 
is still possible below the noise-limit scale / min but for observational 
data, this may not be reliable any more (see Sect. 14.3. lb . 
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5. Application of WWCC to the molecular cloud 
G333 

To test the application of the WWCC to observed maps we 
selected molecular line maps of the giant molecular cloud 
G 333 obtained with the Mopra telescope ( Lo et ali2009h . Over 
20 molecular transitions were mapped and analysed through 
the Principal Component Analysis (PCA) to investigate spa¬ 
tial correlations among different line transitions and classify 
molecules into high- and low-density tracers. The PCA quan¬ 
tified the similarity of the different maps for eight significant 
principal components. Each component can be visualized and 
the components follow a hierarchy of sizes, but the PCA is 
not able to clearly identify individual spatial scales at which 
the similarity of different maps is broken or where gradients 
between the mapped structures occur. To prove the value of 
the WWCC in identifying these scales, we selected three of 
the maps from lLo et al.l ( 2009 ) that show a different behaviour 
in the PCA. The maps of 13 CO 1-0 and C 18 0 1-0 show very 
similar eigenvalues for the first five eigenvectors, indicating a 
very good correlation at most scales. The map of HCN 1-0 in 
contrast only shows the same eigenvalue for the first principal 
component, but a strongly deviating behaviour for basically all 
other components. Here, we will try to quantify the differences 
and commonalities seen by the PCA for the three maps in terms 
of the correlation and displacement between structures in the 
maps as a function of their spatial scales. 

To stay comparable to the results of the PCA, we follow 
the same procedure in producing the analysed total intensity 
maps as performed bv lLo etld . ( 2009 ). The maps were spatially 
smoothed with an Gaussian of FWHM=45" = 0.75 arcmiio 
and the molecular lines of 13 CO and C 18 0 were integrated over 
the velocity range from -70 to -40 km s -1 . A wider velocity 
range from -80 to -30 km s -1 was chosen in case of HCN be¬ 
cause of visible line wings. The maps are shown in Fig. [33] At 
the distance of 3.6 kpc of G333 (iLockmanlll979 ), the resolu¬ 
tion of 0.75 arcmin corresponds to 0.79 pc. 

In the next step, we need to generate significance maps 
(weights) for the intensity maps to account for the uncertainty 
of the information in the individual pixels due to observational 
noise, stem ming from variable sensi tivities and finite integra¬ 
tion times. lOssenkopf et al.1 ( 2008bl) proposed to use a satu¬ 


rating linear function of the inverse noise level as weighting 
function, w = min(l, axS/N) with a being an appropriate scal¬ 
ing factor. The choice of the saturation level, i.e. the S/N ratio 
above which we assume that the distortion of the structural in¬ 
formation by noise is negligible, is somewhat arbitrary. Based 
on our numerical experiments in the previous sections and the 
inspection of obviously “noisy structures” in the maps, we use 
a value of S/N=25 here, where we define the signal level as the 
median of the intensity over the map. 

We estimate the noise at each pixel from the emission-free 
channels in the velocity ranges from -80 to -75 km s -1 and 
-27 to -15 km s -1 for 13 CO and C 18 0. For HCN, we use the 
velocity ranges -80 to -70 km s -1 and -30 to -20 km s -1 . 


4 The MOPRA beam at the frequenc y of the three lines is 36". The 
smoothing was performed by iLo et al.1 (2009) to allow for a uniform 
grid for all lines from the OTF observations and a reduced noise. 


The significance maps are given next to the integrated intensity 
maps in Fig. [33j Due to the high signal-to-noise ratio of the 
13 CO observations, the S/N=25 limit means that we consider 
almost all the measured information in the 13 CO map as reli¬ 
able, with the exception of a few pixels at the map boundaries 
that have a lower integration time coverage. For the C 18 0 and 
HCN maps, larger regions with low coverage producing some 
noisy patterns are weighted down by their lower significance in 
the statistical analysis of the A-variance and the WWCC. The 
significance maps basically represent the coverage of the map 
in the mosaicing observation with multiple coverages of the 
individual subfields visible as square areas in the significance 
maps. Saturated (white) strips and areas in the maps of weights 
indicate a very good coverage, leading to a low noise so that 
all details of the lines can be seen. The combination of the mo¬ 
saicing observations leads to a high significance (low noise) 
in overlapping areas and a lower significance (higher noise) at 
the map boundaries. To test the sensitivity of our results on the 
choice of the saturation limit we varied the limit by a factor 
two in both directions but found no measurable impact on the 
A-variance spectra and only a small impact on the WWCC on 
the scale of a few pixels, well within the error bars. 

Figure [34] shows the A-variance spectra for the integrated 
intensity maps of 13 CO, C 18 0, and HCN. The C 18 0 and HCN 
clouds have a dominant scale size of 7.2 and 10 arcmin (7.5 
and 10.5 pc), while the 13 CO emission has a larger dominant 
structure with If ~ 12 arcmin (12.5 pc). This matches approxi¬ 
mately to the size of the largest structures with more than half 
of the peak intensity level in the individual maps. 

Based on the results from Sect. 14.3. II we can also estimate 
the enhanced scale of these maps using the critical scale l c 
(Eq. (|40 b computed from noise-corrected maps (see Appendix 
A for an example of the noise-correction). The gradient spectra 
of the noise-corrected A-variance are shown in Fig. [35]in com¬ 
parison to the spectrum obtained for a pure fBm with £ = 2.5 
convolved with a Gaussian beam of 0.75 arcmin. 

In the overall shape, we find a clear similarity between the 
spectra from the simulated map and the observed maps. On 
small scales, there is a peak due to the blurring of smaller struc¬ 
tures through the telescope beam with a_filter size of 0.75 ar¬ 
cmin as discussed already bv iBensch et al.1 ( 2001 ). This dom¬ 
inates the gradient spectra below 2 arcmin. A pure Gaussian 
beam should provide a monotonous increase towards small 
scales as seen for the HCN data. The drop of the gradient spec¬ 
tra for the 13 CO and C 18 0 at the two smallest lags may be ei¬ 
ther due to an imperfect noise correction or to a relative surplus 
of structures at or close to the resolution limit. Therefore, we 
cannot reliably quantify the structural properties of the maps at 
scales below 1 arcmin. 

Above the scales affected by the beam, we also find a signif¬ 
icant difference between the convolved fBm and the observed 
maps. Instead of the flat spectrum with a drop at scales ap¬ 
proaching the map size for the fBm, the observed maps show 
a clear peak around 5-6 arcmin (5.5 pc) and a quick drop 
of the gradient towards larger scales, providing the roots that 
match the prominent scales of 7-12 arcmin (7-12.5 pc) dis¬ 
cussed above. The secondary peaks in the gradient spectrum 
happen on critical scales ^6.3 arcmin for 13 CO and HCN and 
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^4.8 arcmin for C 18 0. Using Eq. (l40l) this translates into en¬ 
hancement scales of « 12 arcmin ( 13 CO and HCN, 12.5 pc) 
and 9 arcmin (C 18 0, 9.4 pc), differing from the measured pro¬ 
nounced scales of 12, 10, and 7 arcmin by at most two arcmin- 
utes, proving both methods to be approximately consistent. 

We can compare the scales found by the A-variance analy¬ 
sis with the results from the PC A analysis bvlEo et al. J2009I) . 
The large-scale distribution of the emission in all lines is traced 
by the first principal component (their Fig. 14), showing typical 
structure sizes of 5-10', but the PC A only measures to what de¬ 
gree the three maps match to that component. It does not quan¬ 
tify the match or mismatch in terms of a size scale and it does 
not measure the difference in the dominant structure size for 
the three maps. The clump analysis also performed by Lo et al 


( 2009i) only finds clumps that are much smaller, with sizes < 2'. 

Measurements of (g) and g max of a cloud allow to estimate 
the spectral index by interpolating between £ = 2.5, 3, and 3.5 
(Fig. [27]). In the inertial range, we expect to find (g) = g- 2. To 
evaluate (g) we are restricted to the range not affected by beam 
smearing, i.e. to l > 3 arcmin. We use the plateau before the 
second peaks (Fig. [35]) to estimate the spectral indices of the 
cloud emission structure as £c 13 o = 2.8, £ib co = 2.3, £hcn = 
2.5. 


In Figs. [36] and [37] we show the results of the application 
of the WWCC to the three pairs of maps 13 CO and C 18 0, 13 CO 
and HCN, and C 18 0 and HCN, in terms of the cross-correlation 
coefficient spectra and the spectra of displacement vectors. For 
the pair of 13 CO and C 18 0 maps the cross-correlation coeffi¬ 
cient shows a strong match (r(Z) ~ 0.85 - 0.89) for all scales 
> 2 arcmin (Fig. [36]). In spite of the apparent differences in the 
two maps, this confirms the strong similarity of the emission 
in both molecules, tracing the same conditions of interstellar 
gas. The difference at small scales is partially due to the noise 
and partially due to optical depth effects in 13 CO at the densest 
spots in the map. The fact that the correlation does not reach 
unity at larger scales must be basically due to the insufficient 
excitation of C 18 O in low density gas leading to the apparently 
smaller size of the emitting area and to a lesser degree a result 
of relatively small systematic displacement. The WWCC re¬ 
sult, however proves that this effect is scale-independent (above 
2 arcmin), i.e. we find the same amount of thin gas, dark in 
C 18 0 but bright in 13 CO on all spatial scales in the molecular 
cloud. This is fully in line with the fractal description of molec¬ 
ular clouds. 

When comparing the CO isotopes with HCN, however, we 
find a weaker correlation on all scales and in particular a sig¬ 
nificant drop of the CC coefficient at the scale below « 7 ar¬ 
cmin (7 pc). The weakest correlation is found between C 18 0 
and HCN (r(/ > 2 arcmin) « 0.7). The increase of the CC coef¬ 
ficient at 8 arcmin scale is also visible here but less pronounced. 
This weak correlation contradicts the frequently made assump¬ 
tion that both species are high-density tracers, being sensitive to 
the same dense gas concentrated in small clumps and cores. As 
both molecules have very different critical densities our weak 
correlation at small scales probably reflects different excita¬ 
tion conditions there, i.e. a thermalization of the CO isotopes 
over larger scales (lower densities) compared to HCN. As all 
high-density clumps are embedded in lower density envelopes, 


the correlation is still well above 0.5, but significantly different 
from unity. Possible physical reasons for kink of the CC curve 
at the scale of « 7 arcmin will be discussed in Sect. 15.21 


As the PCA is based on the computation of the cross¬ 
correlation matrix (Eq. (flOb ) we can directly compare our scale 
dependent cross-correlation coefficients (Eq. (fl5l ) , Fig. [36]) with 
the global cross correlation coefficients already computed by 
IFo et al.1 ( 2009 ). They measured coefficients of 0.88, 0.79, and 
0.75 for 13 CO and C 18 0, 13 CO and HCN, and C 18 0 and HCN 
maps, respectively, wheji comparing all pixels of their inte¬ 
grated intensity maps. ILo et all ( 2009 ) did not specify the un¬ 
certainty of the coefficients, so that we repeated their computa¬ 
tion of the correlation coefficients and obtained errors of about 
0.01. When we compare these global CC coefficients with the 
average of the spectrum of the CC coefficient r(/) over the full 
range of spatial scales, the scale-averaged CC coefficients are 
slightly lower (0.82+0.04,0.75+0.05, and 0.69+0.05 for 13 CO 
and C 18 0, 13 CO and HCN, and C 18 0 and HCN maps) than the 
global ones computed without weighting function. The differ¬ 
ence is marginally significant and can be due to two effects: 
the inclusion of the weighting function tends to lower the CC 
coefficients on large scales by 0.02-0.05, and the small CC co¬ 
efficients on the noise-dominated scales < 2 arcmin further re¬ 
duce the average cross correlation coefficient. When restricting 
ourselves to scales without significant noise distribution, / > 2 
arcmin, we obtain numbers closer to the direct correlation co¬ 
efficients measured from PCA (0.85, 0.78 and 0.73). This in¬ 
dicates that the global cross-correlation coefficient traces dif¬ 
ferences in the maps on scales independent from noise quite 
well, but it does, of course, not deal with a variable reliabil¬ 
ity of the data and it cannot distinguish between differences 
on different noise-free scales, such as the varying correlation 
with HCN below and above 7'. The WWCC method is capable 
of recovering the CC coefficients, also measured by the PCA, 
and, in addition, it allows the analysis of the CC coefficients on 
scale-by-scale basis. 


Figure [37] shows the distribution of measured displacement 
vectors as a function of scale for the three pairs of map^. There 
is no displacement between the structures of 13 CO and C 18 0 
clouds (Fig. [37](a)) on small scales (< 10 arcmin), while the 
structure is offset at larger scales with an amplitude of ~ 1 ar¬ 
cmin, where all structures in 13 CO are shifted to the south-west 
relative to C 18 O. If we consider the pairs with HCN, however, 
we find a much larger offset at large scales (Fig. [37] (b) and 
(c)). There is a gradual increase of the offset up to « 3.5 ar¬ 
cmin (HCN relative to 13 CO) and up to 5 arcmin (HCN relative 
to C 18 0), where HCN is offset relative to the CO isotopes in 
south-west direction along the filament. 
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5.1. Identification of correlated structures 

To understand the nature of the structures that are actually re¬ 
sponsible for the measured shift and correlation we go one step 
back and look at the individual filtered maps and compare them 
as for the fBm maps in Sect. 14.2. H Based on the outcome of the 
WWCC, we examine here structures with a scale of 5 arcmin, 
i.e. below the kink in the cross-correlation spectrum at 7 ar¬ 
cmin, and structures at a scale of 18 arcmin, i.e. showing the 
maximum offset between the maps, for the pair of C 18 0 and 
HCN maps. 

Figure [38] (top panels) shows the maps wavelet-filtered on 
the scales of 5 arcmin. The correlation coefficient of the maps 
at this scale is about r(l = 5 ) « 0.7 and the offset between 
them is small r(/ = 5 ) ~ 0.5'. In the figure, we can identify 
two relatively bright cores in C 18 O in the central region of the 
cloud, one of them also matching a peak in HCN, while the 
second one breaks up into two knots in HCN. The peak in the 
filtered HCN map is further offset to the south-west, not having 
any direct counterpart in C 18 0. Multiplication of these maps 
(P(x, l = 5 ); see Eq. (l29lk Fig. [38] bottom left) reveals three 
relatively bright correlated features at the locations of the three 
knots seen in HCN. In the northern core (relative R.A. « -2 ', 
Dec « -2 ) we find only a small shift of the HCN peak rela¬ 
tive to the C 18 O peak. In the southern core we find one feature 
(R.A. « -7\ Dec ~ 7 ) that is prominent both in C 18 0 and 
HCN. The third peak (relative R.A. « 0.3' and Dec « -10 ) 
stems from the global peak in the filtered HCN map that has 
no bright counterpart in the C 18 O map, but falls on a relatively 
faint elongated region there. The matching peaks dominate the 
cross correlation while the southern knot drives the measured 
offset. The WWCC result at small scales combines the contri¬ 
bution from all three features. The two matching knots, and in 
particular the northern core, brightest in C 18 0, are very similar 
in both maps, probably caused by warm and dense core on the 
scale of< 8 arcmin, have similar emission profiles in both tran¬ 
sitions. They create a cross-correlation coefficient of about 0.7 
and do not contribute to any measured offset. In the third fea¬ 
ture, HCN is clearly offset from the C 18 0 emission. It hardly 
contributes to the cross-correlation function, but it produces the 
measured small global offset. 

In Fig. [39] we show the same analysis for the C 18 0 and 
HCN maps wavelet-filtered on the scale of 18 arcmin. The 
WWCC showed a cross-correlation coefficient r(l = 18 ) « 0.7 
level and a displacement of the HCN structure by r x = 2.8 
arcmin and r y = -4.5 arcmin in south-west direction relative 
to C 18 0 (Figs. [36l (dashed line) and[37](c)). The filtered maps 
are dominated by a cometary structure with a relatively sim¬ 
ilar shape, but a clear mutual offset. The simple product map 
(Fig. [39] bottom left panel) reflects this shape but the displace¬ 
ment between the maps adds a negative arc-like distribution in 
south-west direction and leads to a low amplitude of the prod¬ 
uct maps of at maximum 1.8. When accounting for the offset 

5 Note that the sum of the recovered displacement vectors on large 
scales does not exactly add up to zero because of the finite map size. 
The boundary treatment tends to lead to a small underestimate of the 
displacement amplitudes on large scales, as was shown for the dis¬ 
placed Gaussian circular structures (see Fig. m. 


in shifting the HCN map by the measured displacement vector 
from the WWCC before computing the product map (bottom 
right panel), the negative arc-like structure disappears, the sig¬ 
nal becomes more concentrated in the center of the cometary 
structure, and the amplitude of the product map grows to 2.1. 
The weakening of the correlation coefficient on large scales 
(Fig. [36] dashed line) can be attributed to the offset between 
C 18 0 and HCN structures which progressively increases to¬ 
wards large scales (Fig. [37](c)). It becomes clear that at large 
scales, the WWCC is dominated by the southern core also seen 
at small scales. This is systematically offset between HCN and 
C 18 0 along the filament and shows a cometary-shape structure 
at that scale. 


5.2. Interpretation 


The A-variance spectra show the largest difference between the 
13 CO map and the C 18 0 map while the HCN map exhibits an 
intermediate behaviour (Fig. [34]). The C 18 0 map has a shallow 
spectrum with the smallest prominent size, i.e. a relative sur¬ 
plus of small-scale structures relative to the other two maps. 
This can be seen by eye in Fig. [33] All emission structures 
and in particular the peaks are smaller than the corresponding 
structures in the other two tracers^This should be mainly an 
optical depth effect. Following IIL o et al.ll ( 2009 ) most clumps in 
the source are optically thick in 13 CO so that one only sees an 
outer “ 13 CO-photosphere” while C 18 0 traces the full centrally- 
peaked column density structure. HCN 1-0 has an intermediate 
optical depth. Consequently, the optical depth effects produce 
exactly the structure blurring discussed in Sect. 14.31 As a result 
the A-variance spectrum becomes steeper and the pronounced 
structures appear larger than in the underlying column density 
profile. The true structure is probably best traced by the C 18 0 
map. 


In contrast to this effect, we find a much weaker spatial cor¬ 
relation between the HCN map and the other two maps com¬ 
pared to the mutual correlation between the two CO isotopo- 
logue maps. 13 CO and C 18 0 emission line maps are strongly 
correlated above the noise scale. This is explained by the fixed 
relative abundance and similar excitation conditions, with a 
critical density of only 2 x 10 3 cm -3 , that lead to a match of 
the locations of all local features. In that sense, the opacity 
blurring is only a second-order effect. The small offset of the 
largest large structures between 13 CO and C 18 0 measured in 
the WWCC indicates that the optical depth is somewhat higher 
in the south of the C 18 O peaks compared to their northern en¬ 
vironment, so that the corresponding larger 13 CO emission pat¬ 
tern is slightly shifted to the south of the C 18 0 peaks. This 
means that the column density profile around the cores should 
have shallower decrease to the south compared to the north. 


The globally lower correlation coefficient between HCN 
and the CO isotopologues indicates a systematically different 
structure on all scales. We find in particular a clear “kink” to¬ 
wards smaller coefficients at scales below 7 arcmin. At those 
scales we see in Fig.[38]that the emission of the two main cores 
in the cloud is very different for HCN and C 18 O. In the northern 
core (-6',+6') the HCN emission is systematically shifted to the 
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south relative to C 18 0 and in the southern core (0,-10') HCN 
even peaks at a location south of the C 18 0 that is not promi¬ 
nent in C 18 O at all. This systematic difference could indicate a 
chemical transition or a density gradient at the measured scale 
of about 7 pc. The HCN emission peaks in the south of the 
column density peaks traced by C 18 0 could stem from an in¬ 
creased HCN abundance, e.g. due to UV illumination (see e.g. 
Fuente et al. 120051) . or from an increased volume density. As 
HCN 1-0 has a critical density of 10 6 cm -3 , it traces much 
denser gas than the CO isotopologues. The measured scale is 
significantly larger than the size of i ndividual clump s in G 333, 
measured bv lWong et ID <l2008h anri lLo et alJ d2009h to fall be- 
tween about 0.4 and 2 pc, i.e. it must reflect a global process 
that affects the whole cloud. As such a process continued se¬ 
quential star-formation was proposed bv lNguven et al.l ( 2015 ). 
It would be in line with the compression of gas by newly 
formed stars. [Wong et al .l (20080 have observed this process on 
a scale of 9 arcmin around G 333.6-0.2, the northern region in 
our map, but that region contributes only to a small degree to 
the statistics in our whole map. The inclusion of more species 
that are separately sensitive either to different densities or to 
chemical changes produced by UV radiation from young stars 
should allow to distinguish between both effects. 

The systematic offset of the HCN structures relative to 
structures of 13 CO and C 18 0 at large scales (Fig.lJ7l(b) and (c)) 
can be understood from the emission filtered patterns on large 
scales shown in Fig. [39] The structures correlated on scales of 
18 arcmin have a cometary-like structure oriented in the same 
north-east direction. The HCN peak at that scale sits at the head 
of the structure. The measured displacement vector follows the 
cometary structure due to the shift of the individual HCN peaks 
towards the south-east of the corresponding C 18 0 peaks. The 
combination of this structure and the derived displacement vec¬ 
tor indicate a global anisotropy either in the density structure 
or in the chemical structure. A density structure that is consis¬ 
tent with the observations asks for a statistically higher den¬ 
sity in the south-west compared to the north-east. This could 
be a global feature or a structure where every core shows a 
low density tail towards the north-east. This would also explain 
the small shift between 13 CO and C 18 0 due to the higher op¬ 
tical depths associated with higher densities in the south-east. 
A global anisotropy in the abundance structure asks for a pre¬ 
ferred direction of the external UV illumination of the cloud. 
This is not impossible, but there are no independent observa¬ 
tional evidences for such a field from the surrounding stars. 

Altogether we see that the combination of the A-variance 
analysis and the WWCC allows in principle to distinguish the 
structural effects of optical depth, excitation, and variable abun¬ 
dance as a function of spatial scale. For a real discrimination 
between the three effects, however, more data need to be in¬ 
cluded in the analysis, covering well-defined ranges of critical 
densities, chemical formation conditions, and optical depths. 

6. Summary 

We developed the weighted wavelet-based cross correlation 
(WWCC) to study the correlation and displacement between 
structural changes in molecular clouds as a function of scale. 


By inheriting th e prope rties of the A-variance analysis from 
Ossenkopf et alJ (2008a) the WWCC method can use a weight¬ 
ing function for each data point in the map to quantify their 
significance or variable signal-to-noise ratios. The weighting 
function also allows us to deal with irregular map boundaries 
and it permits the WWCC analysis through Fast Fourier trans¬ 
form which considerably speeds up the computation. 

The method can be applied to observations of interstel¬ 
lar clouds taken in different tracers or at different velocity 
ranges. The WWCC measures the correlated structural changes 
between different maps as a function of their spatial scale. 
Characteristic scales where the behaviour of different tracers 
starts to match or deviate from each other can be interpreted 
e.g. as scales of chemical or phase transitions, or driving and 
dissipation scales in interstellar turbulence. 

We performed a series of numerical tests of the WWCC for 
simulated maps with well-defined spatial properties. It reveals 
that 


- For structures with a well-defined size, the spectrum of 
cross-correlation (CC) coefficients strongly depends on the 
displacement between structures and ratio of their promi¬ 
nent scales. The CC spectrum turns negative on small scales 
if the overlapping area between the two structures becomes 
small compared to the not overlapping area. The degree 
and scale range of the anticorrelation depends mainly on 
the displacement and much less on the size ratio. Beyond 
the scale of the mutual displacement, the correlation mono- 
tonically increases. For structures without dominant scales, 
the CC spectrum also traces the mutual displacement: it 
reaches a minimum close to the scale of the offset between 
the structures. In case of noisy data, the amplitude of the 
CC spectrum is lowered at small scales, both for correlated 
and anti-correlated structures. 

- The displacement vectors can be reliably recovered at 
scales dominated by actual structural information, not ob¬ 
servational noise. For overall self-similar structures, inde¬ 
pendent of an enhancement of particular scales, there is al¬ 
ways enough structural information on all scales that the 
displacement vector can be robustly recovered even for rel¬ 
atively low signal to noise levels (S/N > 5). In case of 
individual structures with a fixed size, the lower limit for 
the reliable measurement of the displacement can be es¬ 
timated from the minimum of the A-variance spectrum, 
/low = 1.5 /min + 1.5 pix and the upper scale is estimated 
from / upp = -0.38 / p + 65 pix, where / p is the prominent 
scale of the lar gest s tructure. 

Mac Low & Ossenkop fl ( 2000 ) showed that the A-variance 


can be used to measure the scale of individual prominent 
structures in maps. We extended the A-variance analysis to 
allow for such an identification of enhanced structure scales 
also in maps that have structure on all scales, but with some 
scales being slightly more prominent than the rest. By using 
the logarithmic gradient of the A-variance spectrum we can 
use the scale of the steepest gradient / c to measure the size 
of the enhanced structure k = 0.52/ c as long as the critical 
scale falls above the noise limit l m j n . 
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As the results of the WWCC are sometimes non-intuitive, it 
can be useful to inspect the individual wavelet-filtered maps to 
understand the origin of a varying cross-correlation and mea¬ 
sured displacement vectors in terms of the individual features 
in the maps that drive the statistical parameters. The product 
of the maps filtered (and shifted) on a specific scale of inter¬ 
est allows to identify the locations of correlated structures and 
recognize their shapes. 

We analysed molecular line maps ( 13 CO, C 18 0, and HCN) 
of the G 333 giant molecular cloud. The A-variance spectrum 
and its gradients (Eq. (l40l) ) allow to measure prominent scales 
in the maps at 12 arcmin ( 13 CO), at 7-9 arcmin (C 18 0), and at 
10-12 arcmin (HCN). For the changedn turbulent, self-similar 
structure, covering all spatial scales we can estimate the spec¬ 
tral indices of the power spectrum from the A-variance gradi¬ 
ent spectrum as £ = 2.8 ( 13 CO), £ = 2.3 (C 18 0), and £ = 2.5 
(HCN). This indicates a gradual difference between the three 
maps with most small-scale structure and the smallest promi¬ 
nent features in the C 18 0 data and the dominance of larger 
structures and a relatively larger prominent structure in 13 CO. 
The HCN map shows an intermediate characteristics. All in¬ 
dices are lower than expected from simple Kolmogorov turbu¬ 
lence, but in the typical range observed in molecular line ob¬ 
servations of other interstellar clouds dFalgarone et al12007 ). 

Application of the WWCC to these molecular line maps 
shows that: 

- 13 CO and C 18 0 emission line maps are strongly correlated 
above the noise scale in spite of their overall very different 
scaling behaviour. This is expected from a constant relative 
abundance, similar excitation conditions, and the higher op¬ 
tical depth of 13 CO. A fixed abundance ratio and identi¬ 
cal excitation conditions for both species provide a match 
of all local features of the map, but a perfect match of all 
structures is prevented by the larger 13 CO optical depth. It 
leads to a suppression of the high density peaks relative to 
ther lower-density environments in the intensity map so that 
all emission patches appear somewhat larger in 13 CO com¬ 
pared to C 18 . The measured V offset of the large structures 
between 13 CO and C 18 0 in southern direction indicates that 
this optical depth effect is slightly anisotropic, due to a shal¬ 
lower density gradient to the south. 

- Both the HCN- 13 CO and HCN-C 18 0 correlations decrease 
on scales less than 7 arcmin (or 7 pc), where the maps trace 
different structures. On larger scales, all species trace the 
overall structure of the molecular cloud resulting in tighter 
correlations. A chemical differentiation of HCN relative to 
the CO isotopologues or a very different volume density 
structure could be the origin of those differences on scales 
smaller than 7 pc. The offset of the HCN structure relative 
to structures seen in 13 CO and C 18 0 measures the relative 
shift of the HCN emission to the south along the filament. It 
may be produced by a higher HCN abundance in the south, 
due to an overall anisotropic radiation field, or a somewhat 
higher density in that direction. 
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Appendix A: Recovering the A-variance spectrum 
of an observed map 

To recover the enhancement scale in an efBm-like molecular 
cloud, we need to correct the measured A-variance spectrum 
for the noise contribution of the observations. Here, we demon¬ 
strate the procedure for the G333 C 18 0 1-0 map introduced in 
Sect. l5l (Fig. l33l (b)). In a first step we measure the noise dis¬ 
tribution across the map by integrating off-line channels free 
of emission. To avoid the need for a radiometric correction, we 
integrate the noise over the same width of a velocity range that 
is used to integrate the line emission (Av no ise = Avii ne = 30 
km s -1 ). In the C 18 O data cube this is achieved by integrating 
the noise over the velocity ranges from -80 to -65 km s -1 and 
from -35 to -20 km s -1 . The resulting noise map is shown in 
Fig. IA.ll In case of 13 CO and HCN maps, the observed spec¬ 
tra do not contain enough channels free of emission to select 
the same velocity width for integrating the noise and line emis¬ 
sion. To correct for the smaller integration range of the noise 
channels, Av no i se < Avii ne we have to scale the integrated noise 
intensity radiometrically by VAvii n e/Av no i se . 

In the next step, we compute the A-variance spectra of the 
observed C 18 0 map and the noise map (respective dashed and 
dotted lines in Fig. IA.21) . A pure A-variance spectrum corrected 
for the noise (full line) is then recovered by subtracting the A- 
variances of the observed C 18 0 map and the noise. The critical 
scale can be measured from the gradients in the corrected A- 
variance spectrum (Eq. (l34l)). 
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Fig. 30. WWCC for two efBm maps (shown in (a), maximum- 
filtered on scales of 15 pix and 5 pix) generated with a S/N = 5. 
The second efBm is shifted by T; = 14.1 pix (r ijX = 10 pix 
and Tiy = 10 pix). Part (b) shows the A-variance spectra of 
the first map (full line) and the second map (dotted line) with 
marks of critical scales (/ c = 6.4 pix and l c = 3 pix, full and 
dotted vertical lines, respectively) calculated from their noise- 
corrected spectra, (c) the correlation coefficients and (d) the 
measured displacement vectors as a function of scale. 
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Fig. 31. Correlation coefficient as a function of scale for efBm 
maps filtered with maximum filter of size k\ = 5 pix and 
Zf 2 = 15 pix and offset by 0, 5, 10, 15, and 20 pix (full line, 
dotted, dashed, dash-dotted, and dash-dot-dot-doted lines, re¬ 
spectively). The dashed vertical lines denote the size of the fil¬ 
ters. 
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Fig. 32. Measured mean displacement ((r)) vs. given displace¬ 
ment (ji) for S/N = oo, 25, and 5 (filled circles, squares, and 
filled triangles, respectively). The circles and triangles are dis¬ 
placed by 1 pixel from the identity line (dashed line) for a better 
visibility. The error bars represent l<x variations. 




Fig. 33. Integrated intensity maps of 13 CO (a), C^_ 8 0 (b) and 
HCN (c) emission lines of the G 333 taken by Lo et al.1 ( 2009 ) 
(left panels) and significance weights w for the different points 
in the maps (right panels). White color in the weight map indi¬ 
cates a significance of unity corresponding to a median signal- 
to-noise ratio of 25 or above. For the 13 CO map, only pixels 
close to the map boundary have a lower S/N, i.e. w < 1. 
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Fig. 34. A-variance spectra of 13 CO, C 18 0 and HCN emission 
maps (full, dotted and dashed lines, respectively). Prominent 
scales of the 13 CO, C 18 0 and HCN maps are marked by vertical 
full, dotted, and dashed lines, respectively. 



Fig. 35. Gradients of the noise-corrected A-variance spectra 
( g(l )) of 13 CO, C 18 0 and HCN emission maps (full, dotted 
and dashed lines, respectively). Critical scales of the A-variance 
gradients of 13 CO, C 18 0 and HCN are marked by vertical full, 
dotted, and dashed lines, respectively. The dashed-dotted line 
represents the gradient spectrum of a pure fBm map (£ = 2.5) 
convolved with a Gaussian filter of FWHM=0.75 arcmin. 


Fig. 36. Dependence of the CC coefficient on spatial scale for 
three pairs of emission maps, 13 CO and C 18 0, 13 CO and HCN, 
and C 18 0 and HCN (full, dotted and dashed lines, respec¬ 
tively). 



Fig. 37. Dependence of the displacement vector on spatial scale 
for three pairs of emission maps, 13 CO and C 18 0 (a), 13 CO and 
HCN (b), and C 18 0 and HCN (c). 
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Fig. 38. Top panels. Integrated intensity maps of C 18 0 and 
HCN emission lines wavelet-filtered on scale of l = 5 arcmin 
(left and right panels, respectively). Bottom left. The product 
map generated by multiplying the two wavelet-filtered maps 
(top panels). 
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Fig. 39. Top panels. Integrated intensity maps of C 18 0 and 
HCN emission lines wavelet-filtered on scale of / = 18 arcmin 
(left and right panels, respectively). Bottom left. The product 
map generated by multiplying the two wavelet-filtered maps 
(top panels). Bottom right. The shifted product map of the top 
left and top right maps where the latter one is shifted by mea¬ 
sured offset in the North-East direction (r x = -2.8 arcmin and 
r y =4.5 arcmin). Note the increase of the signal for the shifted 
combination compared to the direct product. 
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Fig. A.l. The noise extracted from the C 18 0 emission line map 
of the GMC G333. 
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Fig.A.2. A-variance of the observed C 18 0 emission line map 
(dashed line) and the noise emission across the C 18 0 map in 
Fig. IA.1I (dotted line). Their difference represents a pure A- 
variance free of the noise effect (full line). The vertical line 
denotes the critical scale l c = 0.66 arcmin at which the pure 
A-variance curve is the steepest. l<x error bars are presented. 




























